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Abstract. The concept of statistical solution of the three-dimensional Navier- 
Stokes equations is addressed, and a new type of statistical solution is given. The 
concept of statistical solution have been introduced as a rigorous mathematical ob- 
ject to formalize the notion of ensemble average in the conventional statistical theory 
of turbulence; a better understanding of this concept is of fundamental importance 
for a rigorous mathematical approach to the theory of turbulence. The concept of 
statistical solution, as introduced by Foias and Prodi in the early 1970's, is that 
of a time-dependent family of Borel probability measures in the phase space of the 
system representing the probability distribution of the velocity field of the flow at 
each given time. The new type of statistical solution is obtained as a projection of 
a measure in the trajectory space, akin to the notion introduced by Vishik and Fur- 
sikov in the late 1970's. This new type of statistical solution has a number of useful 
analytical properties, which we study here and which will be further investigated in 
subsequent works. The current work is devoted only to the time-dependent case. 
The particular case of stationary statistical solutions will be addressed elsewhere. 

I. Introduction 

Our aim in this work is to address the mathematical formulation of the concept of 
statistical solution of the three-dimensional Navier-Stokes equations for incompress- 
ible fluids. Statistical solutions have been introduced as a rigorous mathematical 
object to formalize the notion of ensemble average in the conventional statistical 
theory of turbulence. 

In turbulent flows, physical quantities vary rapidly and erratically in space and 
time but are somewhat well-behaved in a statistical sense, when averaged in some 
form. Averages might be taken over a certain time interval, over a certain region 
in space, and over an ensemble of flows (e.g. a number of experiments in a wind 
tunnel under seemingly the same conditions). It is this latter average which is called 
ensemble average. 
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The conventional theory of turbulence has relied in most part on empirical evidence 
and heuristic arguments [39, 41, 42, 30, 31, 4, 27, 36, 26, 35]. More recently, a number 
of rigorous results have been obtained for mean quantities of three-dimensional flows, 
usually either from time averages of weak solutions (e.g. [29, 5, 8, 9, 13, 17]) or from 
statistical solutions (e.g. [28, 16, 19, 20, 3, 38]). It is therefore our belief that a better 
understanding of statistical solutions are of fundamental importance for a rigorous 
mathematical approach to the theory of turbulence. 

There are two main notions of statistical solutions, one introduced by Foias and 
Prodi [15, 21] and the other by Vishik and Fursikov [45] (see also an earlier related 
mathematical work by Hopf [28]). In the present work we essentially formulate a 
modified definition of Vishik-Fursikov statistical solutions in slightly different form 
than their original definition, and which becomes a particular case of a statistical 
solution in the sense of Foias and Prodi which is more amenable to analysis and has 
a number of useful properties. 

A statistical solution as defined in [15, 21] is a family of Borel measures parametrized 
by the time variable and defined on the phase space of the Navier-Stokes equations, 
representing the probability distribution of the velocity field of the flow at each time 
(Definition 3.2). The definition given in [45], in its turn, is that of a single Borel 
measure on the space of trajectories, representing the probability distribution of the 
space-time velocity field (Remark 3.2). 

The phase space considered here for the Navier-Stokes equations, which we denote 
by H, is the space of square-integrable divergence-free velocity fields with the ap- 
propriate boundary conditions. We consider either no-slip boundary conditions on a 
bounded set Q C K 3 with smooth boundary, or periodic boundary conditions on a 
domain Q = II 3 =1 (0,L i ), Lj > 0, i = 1,2,3. It is assumed that the flow is forced by 
a given external field of possibly time-dependent volume forces with values in H and 
essentially bounded in H with respect to the time variable. In the periodic case, it is 
also assumed that the space averages of the velocity and force fields are zero. 

The concept of Leray-Hopf weak solution is essential to our analysis and refers to a 
weak solution in the classical sense of the Navier-Stokes equations which also satisfies 
a certain energy inequality and which is strongly continuous at the initial time. This 
energy inequality and the strong continuity at the initial time play a crucial role in our 
formulations of the statistical solutions. Leray-Hopf weak solutions are also weakly 
continuous at any given time, a fact that led us quite naturally to make extensive use 
of the weak topology of the space H and of the associated Borel measures, allowing for 
the use of a number of results in measure theory which are, however, quite delicate. 
This idea of using the weak continuity of weak solutions in the study of statistical 
solutions is originally due to Prodi (see [21]) and is exploited here extensively to our 
notion of Vishik-Fursikov statistical solution. 

We first define a Vishik-Fursikov measure, which is more akin to their original 
definition of statistical solution, being indeed a measure in trajectory space. More 
precisely, in our case, a Vishik-Fursikov measure is a Borel probability measure in 
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the space of weakly continuous functions and which is carried by the set of Leray- 
Hopf weak solutions (and with finite mean kinetic energy at each time). Then, our 
definition of Vishik-Fursikov statistical solution is that of the family of projections, 
in time, of a Vishik-Fursikov measure. 

This definition of Vishik-Fursikov measure is to be compared with the original 
definition of statistical solution in the sense of Vishik and Fursikov [45, 46]. Their 
definition is that of a measure on the space of trajectories satisfying a less sharp mean 
energy inequality and such that there exists a measurable subset of the set of weak 
solutions not necessarily of Leray-Hopf type and which carries the measure. 

One noticeable difference is that, in their definition, the measurability of the carrier 
is part of the assumption, while in our approach we prove that the set of Leray-Hopf 
weak solutions is measurable. 

The study of the topological structure of the set of Leray-Hopf weak solutions is 
indeed a delicate part of the work, which is related to the fact that it is not known 
whether the weak solutions are unique or not. Among other similar results, we prove, 
in fact, that the set of Leray-Hopf weak solutions on an interval of the form [0, T] is 
a (^-subset of the set of weakly continuous functions on [0, T] which are Leray-Hopf 
weak solutions on the interval (0, T] (hence not necessarily strongly continuous at the 
initial time), and that this latter set is a a-compact subset of the space of weakly 
continuous functions on H, with the topology of uniform weak convergence (Corollary 
2.1 and Proposition 2.5). 

The bounded subsets of the set of Leray-Hopf weak solutions on an interval of the 
form [0, T] is a Borel set, but it is not necessarily compact, so we actually define a 
Vishik-Fursikov measure as a measure in the sequential closure of the set of Leray- 
Hopf weak solutions which have finite mean kinetic energy at each time and which 
are continuous at the initial time in a mean sense (Definition 3.3). The sequential 
closure of the set of Leray-Hopf weak solutions is the cr-compact set mentioned above, 
for which the bounded sets are compact. Compactness is usually a crucial step in the 
proof of existence of solutions of a number of different types of problems and here it 
manifests itself in the application of the Krein-Milman theorem, in which the proof 
of existence of Vishik-Fursikov measures is based (Theorem 3.1). 

Then, another difference in our definition of Vishik-Fursikov measure is that the 
mean energy inequality is not an assumption but actually follows from the fact that 
the measure is carried by the set of Leray-Hopf weak solutions, which themselves 
satisfy an energy inequality. We prove, in fact, a strengthened form of mean energy 
inequality (Theorem 3.2). 

Then, as we mentioned above, projecting a Vishik-Fursikov measure, as defined 
here, to the phase space, at each time, yields a family of measures which is a statistical 
solution in the sense of Foias and Prodi (Theorems 3.3, 3.4, and 3.5). This particular 
type of statistical solution is termed a Vishik-Fursikov statistical solution (Definition 
3.4), and it is much more amenable to analysis. 
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Then, we prove a regularity result saying that a Vishik-Fursikov measure is actually 
carried by the set of Leray-Hopf weak solutions (Theorem 4.1). 

One natural question then is whether all statistical solutions are Vishik-Fursikov 
statistical solutions. This would be true if the weak solutions of the three-dimensional 
Navier-Stokes equations were known to be unique, but since this is not known, the 
answer to this question is not trivial. We give, however, in the case of families of sta- 
tistical solutions with support uniformly bounded in H, an intrinsic characterization 
of Vishik-Fursikov statistical solutions as limits, at each given time, of convex combi- 
nations of Dirac measures carried by finite collections of Leray-Hopf weak solutions 
(Theorem 4.2). 

In this work we only consider time-dependent statistical solutions. The mathe- 
matical framework in the case of stationary statistical solutions and, in particular, of 
Vishik-Fursikov statistical solutions will be further developed in [25]. In that work, 
Vishik-Fursikov stationary statistical solutions play a major role and are shown to 
have a number of good properties. 

We believe the solid framework presented here will be useful not only for the study 
of stationary statistical flows, such as turbulent flows in statistical equilibrium in 
time, but also in the study of time-dependent flows, such as in the case of decaying 
turbulence. The role of Vishik-Fursikov statistical solutions in the study of turbulent 
flows will be given elsewhere. The dependence of Vishik-Fursikov measures and of 
Vishik-Fursikov statistical solutions on parameters will also be presented in future 
works. 

Although the theory presented here has been developed specifically in the context 
of the Navier-Stokes equations, it can serve as a model to treat other equations which 
have similar properties and the potential pathologies of the Navier-Stokes equations 
(e.g. the Boussinesq equations of thermo-hydraulics) , such as mainly the possibility 
of a lack of uniqueness and the weak continuity in the phase space. 

2. Preliminaries 

2.1. The Navier-Stokes equations and the mathematical setting. In this sec- 
tion we recall some classical results about individual solutions of the three-dimensional 
Navier-Stokes equations, for which the reader is referred to [34, 33, 43, 44, 10, 18]. 

We consider the three-dimensional incompressible Navier-Stokes equations, which 
can be written as 

dw 

— i/Vu + (u • V)u + Vp = f, V ■ u = 0. 

at 

The variable u = (ui,«2> W3) denotes the velocity vector field; the term f represents 
the mass density of volume forces applied to the fluid and is assumed given; the 
parameter v > is the kinematic viscosity; and p is the kinematic pressure. We 
denote the space variable by x = (xi, x 2 , x 3 ) and the time variable by t. 
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We allow two types of boundary conditions: periodic and no-slip. In the periodic 
case we assume the flow is periodic with period Lj in each spatial direction Xi, i = 
1,2,3, and we set Q = IT 3 =1 (0,Lj). In this case, since the equations for averages are 
easy to solve (see e.g. [44]), we also assume that the averages of the flow and of the 
forcing term over Q vanish, i.e. 



u(x,t) dx = 0, / f(x,t) dx = 0. 
u Jn 

In the no-slip case, we consider the flow on a bounded domain Q C M 3 with smooth 
boundary dQ (at least of class C 2 ), and it is assumed that u = on dQ. Other 
boundary conditions such as those for channel flows can be treated similarly. 

In either the periodic or the no-slip case one obtains a functional equation formu- 
lation for the time-dependent velocity field u = u(t) corresponding, at each time t, 
to the function x G Q > u(x, t). For the functions spaces, one starts, in the periodic 
case, with the space of test functions 

V = /u = w| ■ w G C 00 ^ 3 ) 3 ' V ■ w = 0, / n w(x) dx = 0, w is periodic! 
[ n ' with period Li in each direction X{. J ' 

while, in the no-slip case, one considers the test functions 

V={^™ 3 ; V-u = 0}, 

where C^°(Q) denotes the space of infinitely-differentiable real- valued functions with 
compact support in Q. 

In each case the space H is defined as the completion of V under the L 2 (Q) 3 norm. 
The space V is the completion of V under the if 1 (f2) 3 norm. We identify H with its 
dual and consider the dual space V, so that V C H C V', with the injections being 
continuous, and each space dense in the following one. In the two cases we consider, 
we have in fact that V is compactly included in H. 

We denote the inner products in H and V respectively by 

v)l2 = L u(x) ' v(x) dx ' v ^ H1 = I 2 £ ' 7)x~ dx ' 

Jil i=l,2,3 * * 



n 



and the associated norms by |u|^2 = (u, u)^ 2 2 , ||u||#i = ((u, u))^ 2 . 
The duality product between V and V is denoted by 



( u > v ) 



vy, 



and coincides with the L 2 inner product if u e H and v e V. The norm in V is 
given by 

(u, \) V ',V 



|U||y/ = SUp 

vev || v || 
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We first define the Stokes operator as an operator A : V — > V, by duality, through, 
the formula 

(iu,v)v",v = ((u,v)) 5 i, Vu,veV. 
The restriction of the operator A to D(A) = {u G V; Au G if} yields an unbounded 
self-adjoint closed operator A|d(A) : D(A) C H — > H which is positive definite and 
with compact inverse. Hence, it has a countable number of eigenvalues {Aj} jG N, 
counted according to their multiplicity, in increasing order, with each eigenvalue Xj 
associated to an eigenfunction w, . The Galerkin projector onto the space spanned by 
the eigenf unctions associated with the first m eigenvalues is denoted by P m . Since we 
are assuming the domain to have a smooth boundary, we have the characterizations 
D(A) = H 2 (Q) D V and ^4|d(A) = —PlhA, where A is the Laplacian and Plh is the 
Leray-Helmholtz orthogonal projector from L 2 (Q) 3 onto H. 
The following Poincare inequality holds: 

Ai|u|i a < ||u||^, VuG^, (2.1) 

where Ai > is the first eigenvalue of the Stokes operator. 
We also consider the trilinear form 

b(u, v, w) = / (u ■ V)v • w dx, u, v,wGV, 
Jn 

continuously defined on V, which also defines, by duality, a bilinear operator B : 
V x V — > V according to 

{B(u, v), w)v',v = ^( u > v > w ) 5 Vu, v, w G V. 
Then, the functional equation takes the form 

y + i/ylu + B(u,u) = f. (2.2) 

Throughout this work we consider weak solutions and statistical solutions on a 
time interval I C R and for the sake of simplicity we make the following standing 
hypothesis on the forcing term: 

feL°°(I,H). (2.3) 

Given a subset X of H, we also denote by X w this subset endowed with the weak 
topology of H. In particular, H w denotes the space H endowed with its weak topology. 
The closed ball of radius R in H is denoted by Bh{R)- Since if is a separable Hilbert 
space, its weak topology is metrizable on bounded sets, and in particular B H (R) W is 
a completely metrizable metric space. 

With this framework set, we have the following definition of a Leray-Hopf weak 
solution. 

Definition 2.1. A (Leray-Hopf) weak solution on a time interval I C M is defined 
as a function u = u(t) on I with values in H and satisfying the following properties: 

(i) ug ££,(/-, fon^c^n- 
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(ii) du/dt e L\H(I\ V); 

(iii) u e C\ oc (I; H w ) , i.e. u is weakly continuous in H, which means t i— > (u(t), v) 
is continuous from I into M., for every v G H ; 

(iv) u satisfies the functional equation (2.2) in the distribution sense on I, with 
values in V 1 ; 

(v) For almost all t' in I , u satisfies the following energy inequality: 

i|u(t)|i 2 +u£ \\u(s)\\ 2 m ds < i|u(Oli> + j\f(s)Ms))v ds, (2.4) 

for all t in I with t > t' . These times t' are characterized as the points of 
strong continuity from the right for u, and their set is of total measure. 

(vi) /// is closed and bounded on the left, with its left end point denoted by t , then 
the solution is strongly continuous in H at to from the right, i.e. u(t) — > u(t ) 
in H as t — » t$ . 

From now on, for notational simplicity, a weak solution will always mean a Leray- 
Hopf weak solution. 

Notice that condition (ii) is actually a consequence of (i) and (iv). Condition (v) 
on the energy inequality can be interchanged with the assumption that u satisfies the 
following energy inequality in the distribution sense on J: 

^|u(t)|i 2 + H|u(t)||^<(f(t),u(t)) L2 . (2.5) 

The allowed times t' in (v) can also be characterized as the Lebesgue points of the 
function t !->■ |u(t)|| 2 , in the sense that 

lim - [ t+T \n(t)\l 2 dt = \u(t')\l 2 . (2.6) 

r^o+ r J t , 

Since t |u(t)|^ 2 is locally integrable, these Lebesgue points form a set of full 
measure. 

Since u belongs to Lf oc (I; V), condition (v) implies, upon use of the Cauchy-Schwarz 
and Poincare inequalities, 

\u(t)\ 2 L2 + v [ H5)||| 1 d s <|u(0li a + ^||f||i- (t ,, t; H)(*-0 > (2-7) 

Jt' VA \ 

for t' and t as in (v). 

It has been proved in [12] that weak solutions satisfy a strengthened form of the 
energy inequality. More precisely, given a weak solution u on an interval /, and 
denoting by I' the set of full measure in I defined by the points of strong continuity 



1 This is equivalent to 

(u(t), v) i2 = (u(«), v) La + j {(f (r), v) L 2 - K(u(r), v)) ff i - 6(u(r), u(r), v)} dr, 
for every t, s in / and every v in V; see e.g. [43, Ch. 3, Section 1] 
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of u from the right (see the condition (v) in the Definition 2.1), it follows that for 
any absolutely continuous, nonnegative, nondecreasing function ip : [0, oo) —> K with 
ip' essentially bounded, the solution u satisfies 

^(|u(t)|i 2 ) + ^%'(|u( S )|| 2 )l|u( S )||^ ds 

<^(|u(t')|i 2 ) + ^% / (|u( S )|i 2 )(f( S ),u( S )) L2 d S , (2.8) 

for all t' in I' and for all t in I with t > t'. The proof in [12] has been given in the case 
of periodic boundary conditions and for f G L 2 (I;H) with / = [0,T], but it can be 
easily adapted to the case of a bounded domain in M 3 with a smooth boundary and 
with more general forces, on an arbitrary interval /, such as the case f G L°°(I,H) 
that we consider. 

By using an appropriate sequence of test functions in the inequality (2.5) (see [18, 
Appendix II.B.l] for the details or [2, Proposition 7.3] for a different proof), one 
deduces that a weak solution on an arbitrary interval I also satisfies 

Ht)\h < \u(t>)\l 2 e-»^ + ^\\f\\l^ >t;H) (l - e^*-*')) , (2.9) 

for almost all t' in / and all t in / with t' < t. The allowed times t' are again the 
points at which the solution is strongly continuous from the right. 

In the case of a weak solution on an interval / closed and bounded on the left with 
left end point to, condition (vi) implies that the point to is a point of strong continuity 
from the right, hence the estimate (2.9) is also valid for the initial time t' = to- 
Let Ro be given by 

Ro = — r— ||f 1 1 !/*>(/;#)■ (2-10) 

The energy estimate (2.9) implies the following invariance property for any ball of 
radius larger than Rq: If u is a weak solution on [to, oo) and R > Ro, then 

u(t ) G B H (R) u(t) G B H (R), Vt > t . (2.11) 

It is well known that given any initial time to and any initial condition u in H, 
there exists at least one global weak solution on [t , oo) satisfying u(t ) = u . See, 
for instance, [10, 18, 33, 43]. 

2.2. Properties of weak solutions. In this section we address some further prop- 
erties of weak solutions. First, in order to study some topological properties of the 
space of weak solutions (see Section 2.5) we will need the following two results. 

Lemma 2.1 (Pasting Lemma). Let be a weak solution on an interval (ti,t2] and 
u^ 2 ) be a weak solution on an interval [t 2 ,t 3 ), with — oo < ti < t 2 < t 3 < oo and 
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u 



p)(t 2 ) = uW(t 2 ). Then the function 



u(f) 



u 



u 



«(t), ti < t < t : 
(2) (0, h<t<t 



2, 



3, 



(2.12) 



a iuea& solution on the interval I = (t l7 t 3 ). 

Lemma 2.2 (Compactness Lemma). Lei {uj}j £ n ^ e sequence of weak solutions on 
some interval I = (to, ti), — oo < to < ^i < oo, and suppose this sequence is uniformly 
bounded in H. Then, there exists a subsequence {uj'jy and a weak solution u on I 
such that Uj> converges to u in H w uniformly on any compact interval in I. 

The proof of Lemma 2.1 is simple (see [24] for some details). The proof of Lemma 
2.2 is classical and follows from uniform (in j) estimates of the type (i) and (ii) of the 
Definition 2.1 and the Arzela-Ascoli Theorem (see e.g. [10, 33, 43, 44]). 

Another result which will be useful is the following criterion for strong continuity 
from the right for a weak solution. 

Lemma 2.3. Let u be a weak solution on an interval /cK. Let t' £ /. Then, the 
following are equivalent 



r^0+ T J t , 

Proof. Assume (1) holds. Let e n be a sequence of positive numbers with e n — > 0. 
Then, there exists a decreasing sequence of positive times r„ — > such that 



T n Jt' 

Hence, for each n £ N, there exists a set /„ C (t',t' + r n ) of positive Lebesgue measure 
and such that 



Since I n is of positive Lebesgue measure, we can find, for each n £ N, a time t' n £ I n 
which is a point of strong continuity of u from the right. Hence, the following energy 
inequality holds: 




t-)>o+ r J t , 

(2) u is strongly continous from the right at t' ; and 
1 r t '+ T 





n(s)\ 2 L2 - \u(t')\ 2 L2 < e, 



for all s £ /, 
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as n — > oo, which means that 

limsup|u(4)|| 2 < |u(0l| a . 

71— >OC 

This, together with the weak continuity of u at t — t', implies that u(t' n ) converges 
strongly in H to u(i'). Then, passing to the limit in the energy inequality, we find 
that 

i|u(t)& + u £ WutfWli ds < 1 -\u(t')\l 2 + j\f(s)Ms)) L * da, 

for all t G /, t > t'. This implies that t — t' is a point of strong continuity of u from 
the right, which proves (2). 

Assuming now that u(t) converges strongly in H to u(t'), then t h-> |u(t)|? 2 — 
|u(t')|^ 2 is continuous at t — t' and vanishes at t — t', which in turn implies (3). Now 
it is clear that (3) implies (1), which completes the proof. □ 



2.3. Elements of measure theory. In this section, we recall a number of basic 
results in functional analysis in order to help the reader whose background is in a 
more applied such as statistical fluid dynamics. 

The statistical solutions defined in this work are Borel probability measures in 
some appropriate topological spaces. In particular, statistical solutions are defined 
as families of Borel probability measures in the phase space of the system. We also 
consider Borel probability measures in time-dependent function spaces. We shall 
therefore collect a number of basic facts from measure theory which are needed in 
this and other works on this subject. For the results mentioned here, the reader is 
referred to [1, 6, 7, 40, 37, 32]. 

In what follows, we will assume that a topological space is a Hausdorff space, in 
which any pair of distinct singletons can be separated by disjoint open sets. 

We recall that the family of Borel sets in a given topological space X is the smallest 
a-algebra containing the open sets in that topological space. We denote the collection 
of Borel sets in X by B(X). If the topological space is a reflexive Banach space then 
the Borel sets in the strong topology coincide with the Borel sets in the weak topology. 
In particular the Borel sets in H coincide with the Borel sets in H w . In this particular 
case, Borel sets in V are also Borel sets in H (see [18, Appendix IV.A.l, p. 213]). 

A measurable space is a pair (X,Ai) where A is a set and M. is a er-algebra of 
subsets of X called the measurable sets. A measure space is a triple (X,M.,/jl) where 
(X,Ai) is a measurable space and \i is a measure defined on the sets in M.. A 
probability space is a measure space (A, Ai, /i) for which \x is positive and fi(X) = 1. 
A measure is said to be complete if any subset of a null set (a measurable set with 
measure zero) is also measurable. 

Given two measurable spaces (X, M) and (Y,Af) and a function / : X — > Y, the 
function / is said to be measurable, or (Ai, A/")-measurable, if f~ 1 (E) e Ai for all 
E eAf. 
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A Borel measure on a topological space X is a measure /i on X defined on the 
cx-algebra of the Borel sets B(X) of X. When the target space is metrizable (and 
the measure in the target space is the Borel measure) then the pointwise limit of a 
sequence of measurable functions from a measurable space into that metrizable space 
is measurable (see [1, Lemma 4.29]). 

Given a Borel measure fi on a topological space X, the a- algebra B^(X) is defined 
as the smallest cr-algebra containing the Borel sets and the subsets of Borel sets of 
/i-measure zero. One can show that E G B^X) if and only if there exists a Borel 
set E B and a subset E N of a Borel set of /i-measure zero such that E = E B U E N . 
This representation of E G B^(X) may not be unique but the /i-measure of E B is 
independent of the representation, so that we can extend the Borel measure /i to 
a complete Borel measure on B^(X) by defining fi(E) = ^(Eb). Such a measure is 
called the Lebesgue extension of the Borel measure and is still denoted by /i. Moreover, 
we call the elements in B^X) as /> measurable sets. The collection B^X) of /i- 
measurable sets is usually larger than the collection of Borel sets B(X). 

Given two topological spaces X and Y and a continuous function / : X — > Y, it 
follows that / is a Borel map in the sense that / _1 (i?) G B(X) for all E G B(Y). Such 
a continuous function is also (B^(X), 23(F) )-measurable, with respect to the extension 
of a Borel measure /i on X. Notice, however, that given two (Lebesgue extensions 
of) Borel measures /i and v on X and Y, respectively, a continuous function / may 
not be measurable from (X, B^X), /i) to (Y, B V (Y) , v) since f~ l (E) may not belong 
to Bn(X) for all E in B U (Y) (just take / to be the identity to get a contradiction). 

A carrier of a measure is any measurable subset of full measure, i.e. its complement 
is of zero measure. The support of a Borel measure is the smallest closed set of full 
measure. 

Given a topological space X, we denote by C^(X) the space of bounded real- valued 
continuous functions on X, and by C C (X) the space of compactly supported real- 
valued continuous functions on X. 

A locally compact Hausdorff space X is a (Hausdorff) topological space with the 
property that every point has a compact neighborhood. The Kakutani-Riesz Rep- 
resentation theorem [6, 40] asserts that a positive linear functional / defined on the 
space of compactly supported continuous real- valued functions C C (X) on a locally 
compact Hausdorff space X can be uniquely represented by a Borel measure ji on X, 
with 



Given two measurable spaces (X,Ai) and (Y,J\f), a measurable function T : X — > 
Y, and a probability measure /i on (X, M), one can define a probability measure \x? 
on Y by the formula (it{E) = [i(T~ l (E)), for all measurable subsets E of Y. The 
measure \ir is called the measure induced from \i by T, and is sometimes denoted T\x 
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or jjT 1 (see [1, Section 13.12]). It follows that 



/ <p(y) di*r(v) = / <p(t(x)) d/i(x), yip e L 1 ^). 



(2.13) 




(2.14) 



If /i is a regular Borel measure (as defined in (2.15)-(2.18) below) on a locally 
compact Hausdorff space X, then C C (X) is dense in L p (/i), for 1 < p < oo [1, Theorem 
13.9], and (2.13) holds for all tp e C C (X). 

In the case X and Y are locally compact topological spaces, a continuous map 
T : X — > Y induces an operator : C C (Y) — > C C (X) given by Lrpip — ip oT. Then, 
regarding a Borel probability measure \x on X as an element of the dual space C C (X)', 
it is natural to view the induced measure T\i as L^fj,, where L* T : C C (X)' — > C C (Y)' is 
the adjoint of the operator Lt- 

In a metrizable topological space X, the following statements concerning two Borel 
probability measures fi and v are equivalent [1, Theorem 15.1]: 



where V is any dense subset of C^X) 

In a metrizable topological space X we say that a net {fi a } a of Borel probability 
measures on X converges weak-star to a Borel probability measure /i on X if 



ji = v -t=r-fi(G) = v(G) for all open sets G 
<^>jj,(F) = v(F) for all closed sets F 





(See [1, Section 15.1].) This convergence is denoted by 




A Borel measure fi on a topological space X is called regular when 
H(E) = inf{;u(0); O D E, O open in X} = inf jii(O), and 



(2.15) 




(2.16) 



(2.17) 
(2.18) 
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The first relation is called upper regularity and the second one, lower regularity. 2 The 
Lebesgue extension of a regular Borel measure is also a regular measure. 

Now if K is a compact subspace of a locally convex topological vector space V 
and E is the set of extremal points of K then the Krein-Milman theorem [14] asserts 
that the closed convex hull coE of E coincides with that of K. In the case X is a 
compact metric space and we take V = C(X)' and K to be the subset V of positive 
linear functionals with norm one, then the weak-star topology in V is metrizable, and 
K is convex and weakly-star compact and is made of the Borel probability measures 
on X. The extremal points E of K are precisely the Dirac measures on X. The 
Krein-Milman theorem applies and yields that any Borel probability measure on X 
is the weak-star limit of convex combinations of Dirac measures on X. 

A space that plays an important role in measure theory is the Polish space, which is 
a separable and complete metrizable space. An important property is that any finite 
Borel measure on a Polish space is regular in the sense above [1, Theorem 12.7]. This 
property will be particularly exploited in a forthcoming paper [25], which considers 
stationary statistical solutions. 

In our case, all the spaces H, V, V are Polish, and so are the bounded, weakly 
closed subsets of H endowed with the weak topology, such as Bh(R) w , R > 0. The 
space H w , however, is not Polish. 

2.4. Time-dependent function spaces. We define some basic "time-dependent" 
function spaces. In what follows, we consider an arbitrary interval / in M. First, we 
consider the spaces C\ oc (I,H w ) and Ci oc (7, B H (R) W ), with R > 0, endowed with the 
topology of uniform weak convergence on compact intervals in /. With this topology, 
since H w is separable, the space C\ oc (I,H w ) is a separable Hausdorff locally convex 
topological vector space (the proof can be adapted from [1, Lemma 3.99] using that 
the topology is that of uniform convergence on compact subintervals, and using that 
bounded sets in H w are metrizable), and C\ oc (I, Bh(R) w ) is a Polish space. In the 
case / is compact, we write simply C(I,H W ) and C{1 , Bh(R) w )- 

The topology of Ci oc (I, H w ) can be characterized by a basis of neighborhoods of the 
origin, given by 



where J is a compact subset of I (or just a compact interval in /, the generated 
topology is the same) and O w is a (weak) neighborhood of the origin in H w . 
For intervals J C / C M, we define the restriction operator given by 



0(J,O w ) = {v e C loc (l,H w ); v(t) e O w , Vt e J} 



(2.19) 



rij : C\oc{I, -Hw) — > Ci oc (J,H w ) 
u ^ (IIju)(t) = 



u(t), Vt e J. 



(2.20) 



Some authors refer to this notion of lower regularity as tightness, and define lower regularity in 
terms of closed sets. In this sense, any Borel measure on a metrizable space is regular [1, Theorem 
12.5], but not necessarily tight. 
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For the sake of notational simplicity, we do not make explicit the dependence of the 
operator on /. This should be clear in the context. These operators are continuous. 
In case J is closed in J, ITj is also surjective and open, in the sense of taking an open 
set in C\ oc (I, H w ) into an open set in C loc ( J, H w ). 

For each interval / cK and each t G I, we also define the projection operators 

n< : Ci oc (7, H w ) — > H w 

u h-> ll t u = u(i), 

which are also continuous and open, as well as surjective. 

These operators are crucial for passing from trajectory to phase space and are also 
important for the characterization of certain spaces and for proving some topological 
properties for them. 

For example, given u G C\ oc (I,H w ), we have that u is bounded on each compact 
interval J C I. Hence, given any sequence of numbers Rk > 0, Rk — > oo, and any 
sequence of compact intervals J n <5= J, I = U„J n , we can write 

C loc (I,H w ) = f]\JU-^C(J n ,B H (R k ) w ). 

n k 

The bounded sets in C\ oc (I, H w ) can be characterized as the sets B for which there 
exists an increasing sequence of compact intervals J n (1 / with / = U n J n , and an 
increasing sequence of numbers R n > with R n — > oo, such that 

Bcf]U^e(J n ,B H (R n )). 

n 

Any closed bounded set B endowed with the topology inherited from C\ oc (I,H w ) is 
a Polish space, hence C\ oc (I,H w ) can be termed a "quasi-Polish" space (we call a 
quasi-Polish space any topological vector space such that any closed bounded subset 
is Polish). 

The compact subsets of Ci oc (J, H w ) can be characterized as the sets JC for which for 
every compact interval Jc/, the subset Hj/C is equi-bounded with respect to the 
norm of H, i.e. there exists R > such that |u(t)|^2 < R for all u e K, and all t G J, 
and ITj/C is equicontinuous with respect to the uniform structure of C(J, H w ), i.e. for 
any (weakly) open neighborhood O w of the origin in H w , there exists 5 > such that 
u(t 2 ) — u(tx) G O w for every ti,t 2 G J with |t 2 — ti| < <5. 

2.5. Trajectory spaces. We now consider some spaces of weak solutions and study 
their topological properties. In what follows, we consider R > Rq, where i?o is given 
by (2.10), and intervals Id, and denote by 7° the interior of I. 
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{u G C\ 0C (I, H w ); u is a weak solution on J}, 
{u G C\ oc (I, Bh(R)w)', u is a weak solution on /} , 
{u G C\ oc (I, H w ); u is a weak solution on 1°}, 
{u G Ci oc (J, B H (R) W ); u is a weak solution on 1°} 



(2.22) 
(2.23) 
(2.24) 
(2.25) 



endowed with the topology inherited from Ci oc (J, H w ). 

Remark 2.1. It is straighforward to see that Uj C tl\ and Ui{R) C U\{R). When I 
is open on the left, then actually lAi = U\ and Ui(R) = ld\{R). The main difference 
between the spaces Uj and Ltj (and between Uj(R) and U\{R)) is in the case in which 
/ is closed and bounded on the left, for which the solutions in lAj are necessarily 
strongly continuous from the right at the left end point of /, while the solutions in 
Uj are weakly continuous but not necessarily strongly continuous. The space U\ is 
in fact the sequential closure of Uj. In the bounded case, since U\(R) is metrizable 
and complete (see Proposition 2.1), we have that this space is in fact the closure of 



In order to illustrate that the inclusions could be strict when I is closed and bounded 
on the left, or at least that we cannot prove equality, consider the case in which 
/ = [0,T). Let u = u n be a sequence of weak solutions of (2.2) on I = [0,T) with 
initial data u(0) = u 0n . Assume that, as n — > oo, u 0n converges weakly to u in H. 
Then, it is not difficult to deduce (similarly to Lemma 2.2 or see the existence results 
in [10, 18, 33, 43, 44]) that a subsequence of u n converges to a limit u in L 2 (0,T; V) 
weakly and in L°°(0, T; H) weak-star and that u is a weak solution on (0, T) and not 
necessarily on [0, T); that is we are not able to prove (2.4) for t' = 0. The sequence 
{u n } n belongs to Ui but we can only guarantee that the limit point u belongs to 



We start with the following result. 

Lemma 2.4. Let I C K be an arbitrary interval and R > R . Then, Ui(R) and 
U\(R) are not empty. 

Proof. First, since Uj(R) C Uj(R), it suffices to show that Ui{R) is not empty. This 
amounts to showing that there exists a Leray-Hopf weak solution on I with u(t) G 
Bh{R) for all t G /. If / is closed and bounded on the left with left end point t , 
then we take an arbitrary initial condition u G Bh{R), use the well-known existence 
result (see Section 2.1) to obtain a global Leray-Hopf weak solution u on [t , oo), use 
the energy equation and, in particular, the invariance property (2.11) to deduce that 
u(t) G Bh(R) for all t > t and finally restrict u to / to have an element in Ui(R). 

If / is open on the left with left end point to > — oo, then we take a sequence 
of initial times t 0n — > t and take a fixed initial condition u G B H (R) (or even a 



Ut{R). 




□ 
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sequence of initial conditions provided they are all within Bh(R))- Then, we construct 
a sequence of weak solutions u n , each of them in U\t 0njOO ){R) ■ Using Lemma 2.2 and 
a diagonalization process, we find a subsequence converging to a weak solution u on 
(to,oo), with u(t) G Bh{R), for all t > t . Restricting this solution to / yields an 
element m.Ui(R). □ 

Proposition 2.1. Let I cK be an arbitrary interval and let R > Rq. Then, 

(1) The spaces Uj and li\ are separable Hausdorff spaces; 

(2) The space Uj(R) is a separable metrizable space; 

(3) The space U\(R) is a Polish space. 

Proof. Properties (1) and (2) are not difficult to check. The only step which is not 
immediate is the completeness oiU\(R). Since this space is metrizable, it suffices to 
check sequential completeness. The sequential completeness follows then from Lemma 
2.2. □ 

The Leray-Hopf weak solutions belong to Uj, so this is the natural space to consider, 
but the larger space U\ is needed because each space U]{R) is compact as we will 
prove below 3 . Nevertheless, we show that lAi and Uj(R) are at least Borel subsets of 
C ioc (I,H w ) and Cy oc (I, B H (R) W ), respectively. 

Proposition 2.2. Let I C R be an arbitrary interval and let R > R . Then l(\(R) 
is compact in C\ oc (I, Bfj(R) w ) and, hence, it is a compact metric space, and it is 
compactly embedded in L\ oc (I; H). 

Proof. The classical a priori estimates for the weak solutions (see e.g. [10, 18, 33, 43]) 
yield that 

(1) a bounded subset of U\ is bounded in Lf oc (I; V), 



(2) 



|— ; u G bounded subset of W)| is bounded in L^(I; V'). 



Therefore, U\{R) is bounded in Cj, where 

Cj = {u e Ll c (i;V); u'eLtl 3 c (i;V)}, 

endowed with the natural metric. Notice that Cj is a separable Frechet space, hence 
also a Polish space. Moreover, by the Aubin Compactness Theorem, Cj is compactly 
embedded in Lf oc (I; H). Therefore, ti\{R) is relatively compact in Lf oc (I; H). 

By the Arzela-Ascoli Theorem, the sets B = {u G Bi,u' G B2}, where B\ is 
any bounded set in C\ oc (T, H w ) and B 2 is any bounded set in Lf^(I;V), are rela- 
tively compact in Ci oc (J; H w ). From this we infer that U\(R) is relatively compact in 



3 In [18] we said that Ui(R) is complete, which might not be true. Therefore, in some proofs in 
[18], one has to replace Ui(R) with Lij(R), as it is done in the present paper. 
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Since U](R) is closed and included in C\ oc (I, Bh(R) w ), we deduce that Uj(R) is 
compact in C\ oc (I, Bh(R) w ) and, hence, it is a compact metrizable space with re- 
spect to the topology inherited from C loc (7; if w ). It is also compactly embedded in 

LL(i;H). □ 

The next result shows, in particular, that, in case I is an interval open on the left, 
then Uj and U\ are Borel sets in Ci oc (J, H w ), i.e. they are countable intersections 
of countable unions of closed sets in Ci oc (J, H w ). 

Proposition 2.3. Let I C R be an interval open on the left. Then, for any sequence 
{-RfcjfceN of positive numbers with R k > R and Rk — > oo and for any sequence 
{ J n }neN °f compact intervals in I with I = U n J n , we have the characterization 

U I =U\ = f}{jYlfu] n (Rk)- (2-26) 

n k 

In particular, Uj and U\ are Borel sets in C\ oc (I, H w ). 

Proof. Let u 6 Ui = tl\. Since u is weakly continuous, it follows from the Banach- 
Steinhaus theorem that u(£) is bounded in H uniformly for t on any compact interval 
in I. Therefore, for each n, there exists k sufficiently large such that Ilj n u G Uj (Rk)- 
Hence, for each n, u belongs to \J k ILjWj (R k ). This implies (2.26). 

Since each U\(Rk) is compact, and ITj is continuous, the characterization (2.26) 
shows that the spaces are Borel sets in C\ oc (I, H w ). □ 

Proposition 2.4. Let I C R be an interval open on the left. Then, the bounded sub- 
sets ofUi = Uj can be characterized as the sets B for which there exists an increasing 
sequence of compact intervals J n <s / with I = U n J n , and an increasing sequence of 
numbers R n > Rq with R n — > oo, such that 



Moreover, any closed bounded subset of Uj = Uj is compact, and, thus, Uj = U\ is 
quasi- compact, i.e. a space such that every closed bounded subset is compact. 

Proof. The characterization is not difficult to see. The compactness of the closed 
bounded sets follows from the compactness of each set Uj n (R n ) and by using a diag- 
onalization process. □ 

Proposition 2.5. Let I C R be an interval closed and bounded on the left. Then for 
any sequence {Rk}ken of positive numbers with Rk > -Ro o,nd Rk —> oo, we have the 
representation 

U\ = \Juj(R k ). (2.27) 
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In particular, U\ is a-compact in C\ oc (I,H w ), i.e. it is a countable union of compact 
sets in Ci oc (I,H w ). Moreover, any bounded subset oflA\ must be included in Uj(Rk) 
for k sufficiently large. 

Proof. Let u £ U\. Since u is weakly continuous on J, it follows from the Banach- 
Steinhauss theorem that u(t) is bounded in H for t in any compact interval in I 
containing the left end point of /. Then, applying the energy estimate (2.9) starting 
from any "good" point t' in this compact interval, we find that u is bounded on all 
/. Hence, it must belong to Uj(Rk) for k sufficiently large. Therefore, (2.27) holds. 
Since each U\{Rk) is compact, we deduce that U\ is cr-compact. The characterization 
of the bounded sets is easy to see. □ 

Proposition 2.6. Let I C K be an interval closed and bounded on the left and let 
R > R . Then, Uj and Uj(R) are Borel sets in C\ oc (I,H w ). Moreover, for any 
sequence {Rk}k&n of positive numbers with R^ — > oo, we have the characterization 

Uj = UWjOR*). (2-28) 

k 

Furthermore, any bounded subset ofUi must be included in Ui(Rt) for k sufficiently 
large. 

Proof. Let to denote the left end point of I. Denote by T the topology defined earlier 
for C\ oc (I, B H (R) W ) and denote by T the topology in Ci oc (7, B H (R) W ) associated with 
the uniform weak convergence on compact intervals in / and strong convergence at 
to- More precisely, this topology can be characterized by a system of neighborhoods 
of the origin given by 

Af(J,O w ,e) = {v e C loc (I,B H (R) w ); v(t) £ O w , Vt £ J, |v(t )|^ < e, }, 

where J C / is compact, O w is a neighborhood of the origin in B H (R) W , and e > 0. 
This space is metrizable and, in fact, it is a Polish space. 

The topology T is finer than T, i.e. T is contained in T. We claim that the 
corresponding Borel sets generated by these two topologies are the same. Since T 
is finer than T, it suffices to show that every open set in the former topology is a 
Borel set in the latter topology. For that, it suffices to show that every neighborhood 
J\f(J,O w ,e) is a Borel set in C\ oc (I , B H (R) W ). In fact, Af(J,O w ,e) is a Qs set in 
C\ oc (I, Bh(R)w) since it can be written as 

Af(J,O w ,s)= p| (O(J,O w )nO({t },O^)), 

where 0™' £ = {v £ H; |P m v| L 2 < e} and P m is the Galerkin projector defined in 
Section 2.1. 

Consider now Ui(R) with respect to the topology T. As such, it is a closed set 
since the strong convergence at t together with the energy inequality (2.9) with 
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t' = to implies that the limit of a sequence of weak solutions still satisfies the energy 
inequality with t' = t , from which we deduce that the limit solution is strongly 
continuous from the right at t . Note that since the topology T in Ci oc (J, B H (R) W ) is 
also metrizable, it suffices to consider sequences. 

Now, since Ui(R) is closed for T it is Borel with respect to T, hence it is also Borel 
in Ci oc (I,B H (R) w ). 

Now, due to the energy estimate (2.9), any u G Uj is uniformly bounded on /, 
hence it belongs to Ui(Rk) for sufficiently large k. Thus, the characterization (2.28) 
follows and implies that Uj is also Borel. 

The characterization of the bounded sets is not difficult to see. □ 



2.6. Further topological properties for some trajectory spaces. We now show 
an interesting result saying that Uj is a (/^-set within Ll\ (i.e. Uj is a countable union 
of sets which are countable intersections of relatively open sets in Uj). This result is a 
refinement of the result in Lemma 2.6, but the previous one is also presented because 
of its simplicity. We also show, in the case the forcing term is time-independent, that 
Ui is a large set within L(j, in the sense of being also dense in it. First, let us prove 
the following lemma. 

Lemma 2.5. Consider a time interval I which is closed and bounded on the left, with 
left end point t . Then, the function 

*(u,t) = -i-f (|u( S )|£ 2 - |u(i )|| 2 ) ds (2.29) 

is a Borel function in (u, t) G C\ oc (I,H w ) x (/ \ {to})- Moreover, for each fixed 
u G Ci oc (J, H w ), the function t i— >■ \I/(u, t) is continuous in t G / \ {to}- Furthermore, 
the function u i— > liminf t _ ) . t + \l/(u, t) is a Borel function on C\ oc (I, H w ) with 

liminf *(u,t) > 0, for all u G C loc (I,H w ). (2.30) 

Finally, given uGWj, the following statements are equivalent 

(1) liminf^ *(u,t) = 0, 

(2) u G lAi; and 

(3) lim t ^ o+ *(u,t) = 0. 

Proof. For m G N, let P m be the Galerkin projector onto the first m modes of the 
Stokes operator. Consider the functions 
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which are clearly continuous functions on C\ oc (I,H w ) x (I \ {to}). Since u is weakly- 
continuous on /, the function 

t^ \PrnU(t)\ 2 L2 - \P rn u(t )\ 2 L2 

is uniformly bounded on [to,t] and converges pointwise to the function 

<^|u(*)i! 3 -|u(f )l! a> 

as m — > oo. Therefore, by the Lebesgue Dominated Convergence Theorem, the con- 
tinuous real-valued functions \I/ m (u, t) converge pointwise in (u, t) £ Ci oc (I,H w ) x 
(/ \ {to}) to the function \l/(u, t). Thus, \l/(u, t) is a Borel function in (u, t) £ 
C loc (I,H w ) x (J\0o». 

Now for each u £ C\ oc (I, H w ), we note again that 

t ^ \u(t)\ 2 L2 - \u(t )\ 2 L2 

is uniformly bounded on compact intervals of /, so that 

i r* 

tf(u,t) = — — / {\u{s)\ 2 L2 - \u{t )\ 2 L2 ) da, tel,t> t , 

is continuous in £ £ / \ {to}- 

Since u £ C\ oc (I, H w ) is in particular weakly continuous at t , we have that 

|u(t )| < liminf |u(t)|, 

which implies that 

liminf *(u,£) > 0, 

proving (2.30). 

The proof that u h> liminf \l/(u, t) is Borel in u £ Ci oc (I, H w ) uses the fact that 
^ (u, t) is continuous in t £ / \ {to} for each u fixed. Note that 

liminf ^(u, t) = lim inf \l/(u, s)= lim inf \l/(u, s), 

t-»-t+ t-^Oto<s<t n->oo t <s<l/n 

and since \&(u, t) is continuous in t £ 7 \ {to} we have 

liminf Wu, t) = lim inf \l/(u,s). 

t _rt+ n-Kn <s<t, seQ 

Since (t ,t] D Q is countable, the function inf to<s <t ) sS q ^(u, s) is Borel, and hence 
u H- liminf t _ >t + \l/(u, t) is the limit of a sequence of Borel functions, hence it is also 
Borel. 

The equivalence between (1), (2), and (3) follows directly from Lemma 2.3. □ 



Proposition 2.7. Let I C HI be an interval closed and bounded on the left. Then 
Ui(R) is a Qs-set in U\(R), for any R > R . 
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Proof. Let to be the left end point of the interval /. We have from Lemma 2.5 that 

liminf \&(u,t) > 0. 

for all u G Uj(R), and that such u belongs to Ui(R) if and only if 

liminf *(u,i) = 0. 

t->t+ 

Hence, we can write the "bad" set Uj{R) \ Uj(R) as 

Uf(R) = U\{R) \Ui(R) = {u G Uf(R);3 5,e > 0,#(u,t) > e,\/t g (t ,t + 5]} . 
We can also write it as a countable union of the form 

U$(R) = |J Uj ,n (R), 

where 

Up n (R) — {u E tl\{R)] tt(u,t) > l/A^VtG (t ,*o + r/n]}, 
where r > is fixed and is such that [t , to + r ] C /• We claim that 

Uj' n (R) is closed in U\{R). (2.31) 

In fact, since U\(R) is metrizable, it suffices to work with sequences Uj G Uj' n (R) 
converging in this metric (i.e. that of C\ oc (I, B H (R) W )) to some solution u G U\{R). 
Since U\{R) is included in Lf oc (I,H) (in fact, compactly embedded, by Proposition 
2.2), it follows that this convergence also takes place in L 2 (to, to + r; H). Since Uj(to) 
converges weakly to u(t ) in H, then 

|u(£ )|l2 < liminf |u,-(t )|. 

j->oo 

This, together with the strong convergence in L 2 (t ,to + r; H), implies that 
^(u,t) > limsup^(u i ,t) > -, W G (t ,*o + r/n]. 

Hence, u G Uj' n (R), which shows that this set is closed in U\{R). Note that since 
U\(R) is compact in C\ oc (I, Bh(R) w ) (from Proposition 2.2), then so are the sets 
U}' n (R). 

Now, the complement U\{R) \ Uj' n (R) is open in U\{R) and we can write Ui{R) = 
U]{R) \ Uf(R) as a countable intersection of open sets 

uj{r)= n (^(^)\^ ,n (^))- 

Hence, Uj(R) is a &-set inWjCR). □ 
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Since XAi = {J k€N Ui(kRo), in the case I is closed and bounded on the left, the 
following corollary holds. 

Corollary 2.1. Let I C R be an interval closed and bounded on the left. Then Uj is 
a Qsa-set inUj. □ 

In the particular case the forcing term f is time-independent we also have the 
following density result. We do not use this density in this article, but it might be 
useful in the study of stationary statistical solutions. 

Proposition 2.8. Let I C R be an interval closed and bounded on the left and let 
R > Rq. Suppose the forcing term is time-independent, with f £ H. Then, Ui(R) is 
a dense Qs-set in U\(R) and XAj is a dense Gs^-set in Uj. 

Proof. That Ui{R) is a C/5-set in Uj(R) and Uj is a G$ a -set in U\ have been proved 
in Proposition 2.7 and Corollary 2.1. We only need to prove the density. Since any 
u £ Uj belongs to some Uj(R) for some R > R sufficiently large, it suffices to prove 
the result for a given R ^ Rq. 

Let then u £ Ll\{R), with R > Ro. We first show the existence of an extension 
of u to a global solution u £ Mf tooo -\{R) 011 [^o>°°). I n the case I is closed on the 
right, say / = [to, , just paste u with a global weak solution on [T, 00) with the 
initial condition u(T), using Lemma 2.1. In the case I = [to, T), take a sequence 
t n of positive times converging to T from the left, and, for each n, consider a global 
weak solutions on [t n , 00) with initial condition u(t n ). Then, from Lemma 2.1, each 
such global weak solution can be concatenated with u on [to>V] to yield a global 
weak solution on [to, 00). Thanks now to Lemma 2.2 this sequence of concatenated 
solutions have a limit point which is a global weak solution u on [t , 00), which must 
agree with u on /. Note that Lemma 2.2 does not apply to [to, 00) since this interval 
is closed on the left, but since the sequence of solutions agree with u near the origin 
we just need to apply this lemma on, say, (t , 00). 

Let now t' n — > be a sequence of positive times which are points of strong con- 
tinuity from the right for u. Define u n (t) = u(t^ — t + t), for t > t - Since f is 
time-independent, u n is also a weak solution on [t , 00) which belongs to U[ t0>oo )(R). 
Let u n £ Ui(R) be the restriction of u n to the time interval I. From the (uniform on 
bounded intervals) weak continuity of u it follows that u n converges to u in Lij(R), 
which shows that Uj(R) is dense in 1A\(R). □ 

3. Time-dependent statistical solutions 

In this section we address the definition and properties of time-dependent statistical 
solutions and of Vishik-Fursikov measures and Vishik-Fursikov statistical solutions of 
the Navier-Stokes equations, casting a new perspective on the theories given in the 
fundamental works [15, 21, 45, 46] (see also [18]). 
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3.1. Cylindrical test functions. For the definition of statistical solutions one needs 
to consider appropriate test functions in C(H W ). For this purpose we use the following 
definition. 

Definition 3.1. The cylindrical test functions are the Junctionals $ : H — > R of the 

form 

$(u) =0((u,vi),...,(u,v fc )), (3.1) 
where k £ N, <f) is a C 1 real-valued function on R fc with compact support, and 
Vi,...,Vjfc belong to V. For such $, we denote by $' its Frechet derivative in H, 
which has the form 

k 

$'(u) = d j<f>(( u , vi), . . . , (u, v fc ))v fc , 

where dj<j) is the derivative of with respect to its j-th coordinate. 

The significance of the set of cylindrical test functions can be seen by the fact 
that, when restriced to a bounded ball B H (R) W , R > 0, the cylindrical test functions 
are dense in the space C(Bh(R) w )- In fact, since Bh(R) w is a compact and sepa- 
rable metrizable space, and the cylindrical test functions form a subalgebra which 
contains the unit element and separates points in Bh{R) w , if follows from the Stone- 
Weierstrass Theorem [14] that the space of cylindrical test functions restricted to 
B H (R) W is dense in C(B H (R) W ). 

3.2. Definition and existence of time-dependent statistical solutions. Time- 
dependent statistical solutions in the sense of [15, 21] are defined in the following 
way. 

Definition 3.2. For a given interval I C R, a family {fi t }tei of Borel probability 
measures on H is called a statistical solution of the Navier-Stokes equations over I if 
the following conditions hold: 
(i) The function 

/ $(u) d/x t (u) 
Jh 

is measurable on I for every bounded and continuous real-valued function $ 
on H; 



[ii) The function 



belongs to L™ C (I); 
(iii) The function 



belongs to L 1 1 oc (J); 



1 1-> / |u||a d/z t (u) 
Jh 

t^- ||u||#i d)Ut(u) 
Jh 



H 
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(iv) For any cylindrical test function $, the Liouville-type equation 

[ $(u) d/i t (u) = / $(u)d/v(u)+ / / (F(u),$'(u)) nv <Wu)ds (3.2) 

iff JH Jt' JH 

holds for all t', t £ J , where F(u) = f — za4u — -B(u, u), so t/iat 
(F(u), $'(u))^,y = (f , $'(u)) L2 - u((u, &(u)))hi - 6(u, u, $'(u)); 

(v) T/ie strengthened mean energy inequality holds on I, i.e. there exists a set 
I' G I of full measure in I such that for any nonnegative, nondecreasing, 
continuously- differentiable real-valued function ip : [0, oo) — > IR with bounded 
derivative, the inequality 

iKHi,) <Wu) + v f [ ^(MiOH* d^(u)d s 

Jt' JH 

<\( V(|u|i 2 )d^(u)+ f /V(|u|| 2 )(f,u) i2 d// s (u)d S (3.3) 

1 JH Jt' JH 

holds for any t' £ I' and for all t £ I with t > t' ; 

(vi) If I is closed and bounded on the left with left end point to, then the function 

/ ^(|u|| 2 ) d/i t (u) 

JH 

is continuous at t = to from the right, for every function ip as in (v) . 

Remark 3.1. We have the following consequences and remarks concerning the defi- 
nition above. 

(a) Note that (i), (ii), (iii), (iv) imply that the function t — > J H $(u) d/^(u) is 
continuous on / for all $ of the type considered in (iv). In particular, this 
continuity is valid at the initial time if the initial time belongs to the interval 
/, i.e. if / is closed and bounded on the left. 

(b) Note that condition (iii) actually follows from conditions (i), (ii), and (v) with 
ip(r) = r, r > 0. 

(c) If / is closed and bounded on the left, and assuming (i), (ii), (iii), (iv), (v) 
then condition (vi) above is equivalent to assuming that (v) holds for t' — to. 
Indeed, it is immediate to see, using (i), (ii), (iii) and letting t! — > t in (v), 
that (vi) implies (v) with t' = to- 

For the converse, let P m be the Galerkin projector on the first m modes 
of the Stokes operator and let ip £ be a continuously-differentiable approxima- 
tion of ip with compact support and with < ip £ < ip. Then, $ mj£ (u) = 
ip £ (\P m u\ 2 L2 ) satisfies the requirements in (iv), and 

lim / i/j £ (\P m u\ 2 L2 ) d/i t (u) = / ^ E (|P m u|| 2 ) d/^ (u). 

*->*o JH JH 
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Using this relation and the facts that ip £ < ip and |P m u|^2 < |u|^2 in H (so 
that in particular we can apply the Monotone Convergence Theorem [14, 40]) 
we find 

/ ^(MlO aX,(u) = lim / ip £ (\P m u\ 2 L2 ) d/^ (u) 

lim lim / ifj £ (\P m u\ 2 L2 ) d/x t (u) < liminf / ip(\u\ 2 L2 ) d/i t (u). 



On the other hand, using (v) with t' = to and letting t — > t$ , we find 
limsup / ijj(\u\ 2 L2 ) d/i t (u) < / -0(|u|£ 2 ) <W,(u). 



Therefore, the equality holds and (vi) is proved. 

(d) Taking ip{r) = r in condition (vi) we find that 

lim / \u\ 2 L 2 d/i t (u) = / \u\ 2 L 2 d/i to (u). (3.4) 

Then, using (ii), we find that the initial condition necessarily has finite mean 
kinetic energy: 

/ |u|^ 2 d//t (u) < oo. (3.5) 

JH 

(e) Using the arguments in (b) and (c), one can show in fact that J ip(u) d/i^u) 
and J \u\ 2 L2 d/i t (u) are continuous from the right almost everywhere in /, i.e. 
at any point t' allowed in (v). 

(f) A weaker form of the energy inequality (replacing (v)) is sometimes used in 
the definition of statistical solution, namely 

-/ |u||a d// t (u) + v I I ||u|||i d/i,(u)ds 

* JH Jt> JH 

< - I |u|| 2 d// t /(u) + / / (f,u) L 2 d/i s (u)ds (3.6) 

* JH Jt' JH 

for all t £ / and almost all t' £ I with t' < t. For individual solutions, the two 
corresponding energy inequalities, (2.4) and (2.8), are actually equivalent, as 
proved in [12]. It is not known, however, whether the mean versions (3.6) and 
(3.3) are equivalent. The strengthened version (v) is consistent with the usual 
definition of stationary statistical solutions (see [25]). 

The existence of time- dependent statistical solutions in the sense above was first 
proved in [15] via Galerkin approximation (see Theorem 1 on page 254 and Proposition 
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1 on page 291 in [15]). The existence result can be stated in the following way: Let 
to £ K. and let /i be a Borel Probability measure on H satisfying 



Then, there exists a time- dependent statistical solution {p>t}t>t satisfying p to — z^o- 
We present in the next sections a different proof based on an idea given in [18], 
which yields in fact a statistical solution of a particular type, which we term a Vishik- 
Fursikov statistical solution. 

3.3. Definition and existence of Vishik-Fursikov measures. In the approach of 
Vishik and Fursikov, the statistical solutions are obtained through the help of prob- 
ability measures in suitable trajectory spaces. What makes them measures relevant 
to fluid flows is the condition that they be carried by the space of individual weak 
solutions. We also ask them to have finite mean kinetic energy. Inspired by their 
approach we introduce the following definition. 

Definition 3.3. For a given interval I C M, a Vishik-Fursikov measure over I is 
defined as a Borel probability measure p on the space C\ oc (I,H w ) with the following 
properties 

(i) p is carried by U\; 

(ii) We have 



Ju\ 

(iii) // / is closed and bounded on the left, with left end point to, then for any 
nonnegative, nondecreasing continuously- differentiate real-valued function ip : 
[0, oo) — > M with bounded derivative, we have 



Remark 3.2. The Definition 3.3 for a Vishik-Fursikov measure is inspired by the 
definition of time-dependent statistical solution given by Vishik and Fursikov (see 
[46]), but it is actually not the same. Our definition is stricter in the sense that any 
Vishik-Fursikov measure according to Definition 3.3 is a time-dependent statistical 
solution according to Vishik and Fursikov. Considering for simplicity / = [0, T], 
T > 0, their original definition of time-dependent statistical solution on [0, T] is that 
of a Borel probability measure p on the space Z = L 2 (0, T; H) n C([0, T]; V~ s ), for a 
given s > 2, where V a = D(A a ^ 2 ), asl, and satisfying the following conditions: 

(i) p is carried by the space 



(ii) p is carried by a Borel subset W of Z which is closed in C, and consists of 
weak solutions not necessarily of Leray-Hopf type; 






£ = { u e L 2 (0, T; V) n L°°(0, T; H)\ u' £ L°°(0, T; V~ s ); 
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(iii) p satisfies a weaker form of the mean energy inequality, namely 

J {\u(t)\l 2 + ||u||!} dp(u) < C Qf |u(0)|| 2 dp(u) + 1 

for all t G [0,T], for some constant C > which does not depend on t, and 
where ||u||^ is the natural norm for L. 

It has been proved by Vishik and Fursikov (see [46, Chapter 4]), using Galerkin 
approximations, that for any T > and for every Borel probability po on H with 
finite mean kinetic energy, i.e. 

u\ 2 L 2 dp (u) < oo, 

there exists a time-dependent statistical solution p over the time interval [0, T] with 
initial probability il p = po- This notion of statistical solution will not be used in 
the sequel. 

Remark 3.3. As mentioned in the previous section, the original proof of existence 
of statistical solutions in the sense of Definition 3.2 was also based on Galerkin ap- 
proximations [15]. In [18], a different proof was given, which form the basis of our 
present work. The proof in [18] was based on the Krein-Milman Theorem, in which 
the idea was essentially to prove the existence of a Vishik- Fursikov measure in a sense 
somewhat weaker than Definition 3.3 and then prove that the projections of such a 
measure at each time t yield a family of measures which is a statistical solution. In 
that proof, however, we overlooked the fact that Ui(R) may not be compact, and in 
particular the energy inequality for to at the left end point of the time interval does 
not follow in a trivial way. Using instead the compact space U}(R), the proof of exis- 
tence of a Vishik- Fursikov measure follows very much the steps in [18]; nevertheless, 
we include below this proof for the sake of completeness. The only delicate step and 
which was not proved in [18] is the continuity of the moments of the measure at the 
initial time (expressed by condition (iii) in the Definition 3.3), and for that we need 
to take special care. We also remark on another subtle difference between the result 
in [18], in which the specific construction yielded a family of measures satisfying a 
mean energy inequality, while in our case we prove that any Vishik-Fursikov measure 
satisfies a (strengthened) mean energy inequality, not only the one constructed in the 
proof of existence. 

Theorem 3.1. Let ^GR and let p be a Borel probability measure on H with finite 
mean kinetic energy, i.e. 

vl\\ 2 d/i (u) < oo. 

Then, there exists a Vishik-Fursikov measure p over the time interval I = [t , oo) 
such that U to p = p . 



i- 



f 
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Proof. Let us first consider the case in which /x is carried by a ball Bh(R) with 
R>R 0l i.e. fi (B H (R)) = l. 

Consider the measures P m f^o given by P m fi (E) = /i (P~ 1 -E'), for all Borel sets E 
in H, where P m , m G N, are the Galerkin projectors. By the Lebesgue Dominated 
Convergence Theorem, we have 

/ tp(u) dP m /2 (u) = / <p(P m u) dyu (u) ->■ / y?(u) d/i (u), Vy? G C(B H (R) W ), 

JH JH JH 

(3.7) 

which means that 

P m ^o Akb as m — t- oo. 

Since P m BH(R) w is compact and separable and P m ^o is carried by P m BH(R) w , it 
follows by the Krein-Milman Theorem that, for each m G N, there exist J m,n G N, 
#™' n G R, and u^ n G P m B H (R) such that 

3=1 

and 

/V = 2^°i °<? "n*A*o, as m -> oo. (3.8) 

Since for any nonnegative, nondecreasing, continuously-differentiable real- valued func- 
tion ip : [0, oo) — y M. with bounded derivative, the function defined by 

<p(u) = tp(\P m u\l 2 ) 

belongs to C(Bh(R) w ) for each m G N, we have in particular that 

lim / if)(\P m u\ 2 L2 ) d/i™' n (u) = / ^(|P m u|^ 2 ) dP m fi {u) = / ip(\P m u\ 2 L2 ) d/i (u). 
n ^°° Jh Jh Jh 

Since /i™' n is carried by P m BH(R), we have |P m u| L 2 = |u| L 2 for /i™' n -almost every u, 
so that 

lim / ^(|u||a) d/i[T' n ( u ) = lim / ^(\PmU\h) d/C n ( u ) 

n— ^oo J ^ n— too J 

= [ ^(\P m u\h) dMu) < / d /"o(u), Vm G N. (3.9) 

Since Bh(R)w is compact and metrizable, the space C(Bh(R)w) is separable, and 
there exists a countable dense set {</^}^eN in C(B H (R) W ). Let also {^}^ be a count- 
able dense set, with respect to the uniform topology of C([0, oo),R), in the space 
of nonnegative, nondecreasing, continuously-differentiable real-valued functions with 
bounded derivative. 
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For each p £ N, choose, thanks to (3.7), an index m p £ N such that 

1 



<p £ (u) dP mp /i (u) 



H 



[u) d/i (u) 



n 



< 



2p 



Then, using (3.8) and (3.9), choose n p such that 



H 



(fi(u) d/x p ' p (u) - / ip e (u) dP mp no(u 



H 



< 



2p 



and 



H 



M\"\h)dVo P ' np (u)< I M\*\h) <H(u) + -, 



H 



1 



l,...,p. 



1, ...,p, 



■P- 



(3.10) 



(3.11) 



(3.12) 



Set Hq = (j,™ p ' np . Using (3.10) and (3.11), it follows that for each £ £ N and p > 



<p t (u) d/xg(u) 



H 



(u) d/i (u) 



H 



1 1 1 

~~ 2p 2p p 



(3.13) 



Since {<pe}e is dense in C(Bh(R) w ) an d the measures are probabilities, we find that 
/ <p(u) d/xg(u) -> / y>(u) d/i (u), as p -» oo, Vy? £ C(B H (R) W ), (3.14) 



(3.15) 



which means that 

/!q — 51 /io, as p — > oo. 
Moreover, it follows from (3.13) and the density of that 

limsup / VKMi 2 ) d/i£(u) < / ^(|u|| 2 ) d/i (u), (3.16) 

p^oo Jh JH 

for all nonnegative, nondecreasing, continuously- different iable real-valued functions 
ip : [0, oo) — > R with bounded derivative. We now use properties (3.15), and (3.16), 
along with the fact that p v Q is a convex combination of Dirac deltas (given by (3.8)), 
to construct a Vishik-Fursikov measure. 

For each p £ N, let u p be a weak solution on [i , oo) with initial condition u^(i ) = 
u™J'™ p . Since u^J'" p is bounded in if by R, with > Rq, it follows from the energy 
estimate (2.9) that u p (t) £ B H (R) for aU * > t , p £ N, j = 1, . . . , J m ^ n p. 

Consider now the measures 

jmp.np 

p p= 0™ p > np 5 u P, 

3=1 

which are Borel probability measures on Ci oc ([t , oo), H w ) carried by lil ^(R)- Ac- 
cording to Proposition 2.2, the space i/L ^(R) is a compact metric space. Therefore, 
there exists p £ lih ^ (R) which is the weak-star limit of the sequence (or a subse- 
quence if necessary) of measures {p p } p , i.e. 

(f — ^ p, as p — > oo. (3-17) 
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The measure p is our candidate for the desired Vishik-Fursikov measure. We already 
have that p is carried by tlh ^(R) C Uf tQ so that condition (i) of the Definition 3.3 
of a Vishik-Fursikov measure is satisfied. Since (u, t) i— > \u(t)\ 2 L2 is a Borel function 
on Cioc([*o, oo), H w ) x [t , oo) and 

/ |u(t)|| 2 dp(u) = f \u(t)\l 2 dp(u) < R, 

Jut , Jul JR) 

we see that condition (ii) of the Definition 3.3 is also satisfied. 

Before proving condition (iii), let us show that the initial condition IT fo p = p is 
satisfied. Indeed, we have from (3.15) that Il to p p = p p — ^ p on B H (R) W . On the 
other hand, since Il fo is weakly continuous, we have from (3.17) that IL p p rL p. 
Thus, 

n to p = p - 

We now prove condition (iii). Let ip : [0, oo) — > K. be a nonnegative, nondecreasing, 
continuously-differentiable real-valued function with bounded derivative. Since any 
u G Uh QO )(-R) is weakly continuous at to, it follows from the Fatou Lemma that 

/ V(|u(to)|£ 2 ) dp(u) < liminf / V(|u(t)|£ 2 ) dp(u). (3.18) 

[*0'°°) [to>°°J v ' 

Now, since <p(u) = VKI-^mU^)^) is continuous on ul ooQ JR) for each m G N, we 
have that 



/ 



^(|P m u(0|i 2 ) dp(u) = lim / ^(|P m ii(t)|| 2 ) dpP(u) 

p ^°°Jul JR) 

< liminf / V(|u(t)|| 2 ) dp p (u). 
Using then the Monotone Convergence Theorem we find that 

/ V(|u(t)|| 2 )dp(u)= lim / V(l^u(t)|i 2 )dp(u) 

< liminf / ^(|u(t)|| 2 ) dp p (u). 

P^°° Jut JR) 



Using the fact that pP is actually carried by U[ t0tOO )(R) (since it is a convex combi- 
nation of Dirac deltas carried by weak solutions defined on [to, oo)), the strengthened 
energy inequality (2.8) holds p^-almost everywhere and we find the uniform estimate 



V>(Kt)|| 2 ) < ^(|u(to)i| 2 ) + -^||f |li- (t0fO o^(suP^(r))(t - to) 
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for p p -almost every u. Thus, 



^(|u(t)|i 2 )d^(u)< / V(|u(t )|| 2 )d^(u) 



^\nl^to,oo-,H){^'{r)){t-t ). (3.19) 



Therefore, 



limsup / ^(I u (*)li0 dp(u) 

< lim sup lim inf / tp(\u(t)\ 2 L2 ) dp p (u) 



[t ,oc) V 



< limsup liminf / ^(|u(t )|^ 2 ) dp p (u) 

V H^ R ) (3.20) 

+ -T- II f II i»(t ,oo;«) SU P (* - *0) ) 

= liminf / ^(|u(t )li 2 )d/(u) 



[t ,oo) 



< I V(|u(t )li 2 ) dp(u). 

Putting (3.18) and (3.20) together proves condition (hi). This completes the proof 
in the case in which p is carried by a bounded set in H . Note also that taking ip 
identically 1 in (3.19) and considering the Galerkin projetor P m we have 

/ \P m n(t)\ 2 L2 dpP(u) < I |u(t )li 2 d/(u) + ^-||f ||£. (ld>oo;fl) (f - to). 

[*o>°°) [*o>°°) v ' 

Then, passing to the limit as p — > oo using (3.16) and the fact that u i-» |P m u(t)|| 2 
is continuous, we find 

P m u(t)|| 2 dp(u) < / |u(t )|| 2 dp(u) + J-\\f\\l (t-t ). 

Taking the limit asm oo and using the Monotone Convergence Theorem we obtain 
/ |u(t)|| 2 dp(u)< f |u(t )|| 2 dp(u) 

"f|li»(to,oo;ff)(*-*o), Vt>t - (3.21) 



/ 
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Consider now the case in which po is not carried by any bounded set in H. Then, 
there exists an increasing sequence {R k } k , Rk > Ro, with R k — > oo and such that 
p (A k ) > for every k, where A 1 = B H (R 1 ) and A k = B H (R k )\B H (R k _ 1 ), for k > 2. 
We decompose /i according to 



where p, k is given by Pq(E) = p>o(E D A^), for every Borel set E in H, and fceN. 

Now, since by construction fi(A k ) > 0, we normalize each /!q to a probability 
measure by fio(E) = [1q(E) / )j, (A k ) , for every Borel set E in if. With the procedure 
above, we construct Vishik-Fursikov measures p k on U\(R k ) with Ht p k = P>o(E) and 
satisfying (3.21). Set 



which proves condition (ii). 

It remains to prove condition (iii). Let ip be as in (iii) and assume, without loss of 
generality, that "0(0) = 0. We write 



and note that for a given K £ N, since each p k is a Vishik-Fursikov measure, 




P = ^2^o(A k )p k , 



which is clearly a Borel probability measure. Since each p k is carried by U^ tQ ^{Rk), 
then p is carried by lih which proves condition (i). 
From (3.21) for p k , we deduce that 






K 




K 
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The remaining terms we estimate using (3.21) for p k and the assumption that ip' is 
bounded: 



oo „ 

MAO / ^(|u(t)|i 2 )dp fc (u) 

=K+l JU [t ,oo) 

< jr /i (A)(sup^(r)) ( / |u(t )|| 2 dp fc (u) 
i 1 rem. \ Ju} 



k=K+l V W [t ,co) 



+ ^ll f lli°°(to,oo; J H')(* _ *o) 



= jr (sup^'(r)) ( [ |u(t )|^d/xg(u)+/x (A fc )^||f||icc (t0iOO . fl) (t-t ) 

= sup^'(r) ( / |u(t )lia d/ig(u)+ 

reM \Jh\b h (r k ) 

fj, {H \ B H (R K ))^-\\f\\ 2 L oo(t 0tOO . H ^{t - to^j . 

Since the right hand side above goes to zero uniformly in t, as K goes to infinity, and 
taking (3.23) into consideration we see that 



^(|u(t)|£ 2 )dp(u)^ / V(|u(*o)liOdp(u) 

^ Jul > 

,oo) [tO" 00 ) 

as £ — > to > which proves condition (iii). □ 



3.4. Mean energy inequality for Vishik-Fursikov measures. Notice that in the 
Definition 3.3 of a Vishik-Fursikov measure, there is no explicit condition for some 
sort of mean energy inequality. It turns out that this is hidden in the hypothesis 
that the measure is concentrated on weak solutions, for which the individual energy 
inequality holds. In this direction, we have the following result. 

Proposition 3.1. Let I C K be an arbitrary interval. Let p be a Borel probability 
measure on C\ oc (I, H w ) and suppose p is carried by U\, with 

! |u(t)|| 2 dp(u)<oo, (3.24) 

for almost all t E I. Then, 

/>(0l£ 2 dp(u)eL£ c (J), (3.25) 
Ju\ 

[ \\n(t)\\ 2 Hl dp(u)ELUl), (3.26) 
Ju\ 



34 



C. FOIAS, R. ROSA, AND R. TEMAM 



and the strengthened mean energy inequality holds on I, i.e. there exists a set I 1 of full 
measure in I such that for any nonnegative, nondecreasing, continuously- differentiable 
real-valued function ip : [0, oo) — > M. with bounded derivative, we have, for all t' £ V , 
and all t £ I with t > t' , 

J ui fydmh) + ^V'(iu( S )ii 2 )iiu(s)n^ d,} d P (u) 

<^|^(|u(0li 2 ) + ^V'(|u( S )|i 2 )(f( S ),u( S )) i2 d S } dp(u). (3.27) 

Proof. Let J C / be a closed and bounded interval and recall the restriction operator 
LTj, which is continuous from C\ oc {I, H w ) to C(J, H w ). We know that IijU\ C Uj. We 
argue that it suffices to show the result for the measure pj = Hjp, i.e. that 




< J ui {lHWt')\l>) + ^V'(|u( S )|i 2 )(f( S ),u( S )) i2 d S } d Pj (u) (3.28) 

holds for almost all t' £ J and all t £ J, t > t', , where ^ : [0, oo) — > 1R is an 
arbitrary nonnegative, nondecreasing, continuously-differentiable real- valued function 
with bounded derivative. 

In fact, since Ti'/U) D W} and p is carried by (3.28) implies the corresponding 
result for p and Li\ but still for almost all t' £ J and allt £ J with t >t' . But since 
J C / is an arbitrary closed and bounded interval in J, the inequality extends to 
almost all t' £ / and all t £ / with t > t', which is what we want to prove. Hence, we 
now only need to prove (3.28) 

Since each u £ Uj is weakly continuous and J is compact it follows that u is 
uniformly bounded in H on J. Hence, we can write 

fceN 

Therefore, we first prove the strengthened mean energy inequality in each Uj(R) with 
R > Ro, i-e. that there exists J' of full measure on J such that 

/, \ln<t)\ 2 L ,)+u /VCK^HU* da} d Pj (u) 

< f {^(|u(t')|| 2 )+ /V'(|u( S )H 2 )(f( S ),u( S )) L2 d S } d Pj (u) 

JU S j(R) l Z Jt' J 

holds for all t' £ J' and all t £ J, t > t', and for all ^ as in the statement of the 
proposition. 



TIME-DEPENDENT STATISTICAL SOLUTIONS OF THE 3D NSE 35 

We have seen in Proposition 2.2 that Uj(R) is a compact metric space (with the 
topology inherited from C(J, H w )) and is compactly embedded in Lf oc ( J, H). 

If p(Uj(R)) = the mean energy inequality is trivially satisfied, so we need only 
consider the case in which < p(Uj(R)) < 1, which necessarily happens for R > R 
sufficiently large. 

Let pj be the restriction of the measure p,j to Uj(R) normalized to a probability 
measure (i.e. pj{E) = pj{E C\U](R)) / pj{U]{R)) for any Borel set E in U]{R)). 
Since Uj(R) is compact and separable we can apply the Krein-Milman Theorem to 
approximate fij by finite convex combinations of Dirac measures concentrated on 
weak solutions Uj e Uj(R), j = 1, . . . , J(n), i.e. 

Jin) 

Pn = ^2 V n) P J " as n ~^ 00 ' 

3=1 

with < 6^ < 1, Y^ij=i @j = 1) where the convergence is in the weak-star sense: 

/ V?( u ) dp„(u) / <p(u) dpj(u), 

for all y> e C(U](R)). 

Since each uj™'' is a Leray-Hopf weak solution on J, we have the individual strength- 
ened energy inequality (see (2.8)) 

< ^(luf^Oli^ + ^V'Clu^llaJCf^^f^))^ ds, 

for almost all t' E J and all £ £ J, t > i! . The set jj™^ of allowed times £' above 
depends on the solution (and not on i/j), but since we have a countable family of 

~ (n) 

solutions, the intersection J = C\j, n Jj of the allowed times is still of full measure in 
J, hence the energy inequality above holds for almost all t' in J, independently of 
j, hgN. Thus, considering the convex combination of the solutions we write 




for all t' e J and all t G J with t > t'. 
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Integrate the relation (3.29) in t' to find 



T f \U(\u(t)\j>) + U /V'du^HOIIu^lllxd,} dp B (ll)df 
Jh Ju S j(R) l Z Jt' J 

^ / <2 A {^(l u (Ol! 2 )+ /V / (|u( S )|i 2 )(f(^) ! u( S )) L2 ds) dp n (u)df, (3.30) 



for all ti,t2,t G J, ti < t2 < £. 

Let ti, <2 € J, ti < ^2, and consider the function 



Ju 



which is a well-defined real- valued Borel function on Uj(R), as the limit of the contin- 
uous functions <Pm( u ) = <p(P m ii). Since U\{R) is compactly embedded in Lf oc (J, H) 
(Proposition 2.2), it also follows that <p(u) is a continuous function in Uj(R). Hence, 
since p n converges weak-star to pj we have, using the Fubini Theorem, that 



/ 2 / ^(|u(t)|i 2 ) dp n (u) dt' = [ <p(u) d Pn (u) 

./ti Jul(R) Ju'jiR) 

-> / ^(u) = r / dp-xu) df. 



Consider now the functions 



<p m (u) = J ^(|u( S )|i 2 )||P m u( S )||^ ds. 



Since Uj(R) is compactly embedded in Lf oc (J, H), it follows that <p(u) is a continuous 
function in U\(R). Then, using the Monotone Convergence Theorem and the fact 
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that if)' > 0, we deduce that 

^(|u( S )|i 2 )||u( S )||^ dsdpj(u) 

hm / f ^(|u( S )|| 2 )||P m u( S )||^ ds dpj(u) 



u](R) Jt' 



m— >oo 



lim / (fm(u) dpj(u) 

IU*{R) 



>U}(R) Jt 

for any t, t' G J, t' < t. Similarly. 



= lim lim / y? m (u) dp n (u) 

= lim lim / f iP'(\u(s)\ 2 L2 )\\P m u(s)\\ 2 H1 dsdp n (u) 
< liminf / f ^{\u{s)\l 2 )\\n{s)f m ds dp n (u), 



\ f i)(\u(t)\ 2 L2 ) dpj(u) = lim ~ f i/)(\P m u(t)\ 2 L2 ) dpj(u) 
= lim Urn \ f i/){\P m n{t)\ 2 L2 ) dp n {u) < liminf J / ^(|u(t)|£ 2 ) d/o„(u) 



On the other hand, since the function 



/ il)'(\u(s)\l 2 )(f(s),u(s)) L *ds 
Jt' 

is continuous on Uj(R), it follows that 

/ [ ij)'(\n(s)\l 2 )(f(s),u(s)) L .dsdpj(n) 

JWj{R) Jt' 

= hm / / ^(|u(s)|i 2 )(f(s),u( S )) i2 d S dp re (u). 

n ^°° Ju]{R) Jt' 

Thus, passing to the limit in the mean energy equation for p n we find that 
n 'MHt)\ 2 L2 ) + u f tP'(\u(s)\ 2 l2 )\\u(s)\\ 2 h1 ds) dpj(u) dt> 




'h JU S j(R) I 2 



r / f^(|u(f)l£0 + /V / (|u( S )|i 2 )(f( S ),u( S )) L2 d S } dpj(u)dt, (3.31) 

./ti </«*(/?) l z Jt' J 



< 

for all ti, t 2 ,t G J with ti < t 2 < t 
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The aim now is to divide (3.31) by (t 2 — ti) and take the limit as t 2 —> tf. The 
integrand (with respect to if) of the first term on the left hand side of (3.31) is 
independent of t', while the integrand of the second term in the left hand side and 
that of the second term in the right hand side are continuous in if and present no 
further difficulties either. The only delicate term is the first one on the right hand 
side. 

For that, we consider a countable dense set {ipe}e, with respect to the uniform 
topology of C([0, oo), R), in the space of nonnegative, nondecreasing, continuously- 
different iable real-valued function with bounded derivative. For each £, the function 



Ju](R) 1 



is integrable on J, and hence, by the Lebesgue Differentiation Theorem, there exists 
a set Jg of full measure in J such that for all t\ G Je, 



l —r I' I lMHt)\ 2 L2 ) dpj(u) f lM\<ti)\h) d PA u )> as *2 ti 
t2 - H Ju JuUr) 1 JuUr) a 



■•I 2 

'ti JU" T {R) * JU"j(R) 



Setting J' = f] £ Ji, we have that J' is of full measure and the convergence above holds 
for all ti G J and all £ G N. Thus, dividing the expression (3.31) by t 2 — h, with 
ip = ip£ for each £, and passing to the limit as t 2 — > , we find that 



1 



-MWt)\h) + " jf ^(|u( S )|i 2 )l|u( S )||^d S | dpj(u) 
< f \ 1 -MHti)\h)+ t ^KWs)\lms)Ms))L* ds) dpj(u), (3.32) 

Ju](R) V Z Jt! J 

for all ti G J', alH G J with t\ < t, and all £ G N. Using the density of {ipe}i<=N and 
renaming ti as t', we find 

f \U(\<t)\h) + u t tP'(\u(s)\ 2 l2 )\\u(s)\\ 2 h1 ds) dpj(u) 

Ju S j(R) l Z Jt' ) 

^(|u(t / )| 2 L 2 ) + j("% , (|u( S )l! 2 )(f( S ),u( S )) i2 d S } dpj{nl (3.33) 



ul(R) 12 

for all t' G J', all t G J with ^ < t, and all '0. This proves the strengthened mean 
energy inequality in Uj(R) 
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Applying now the Cauchy-Schwarz, Poincare, and Young inequalities to the forcing 
term in the mean energy inequality (3.33) with ip(r) = r, we find that 

/ \\u(t)\h + u t \\u(s)f m dsl dpj(u) 

< f |u(t')|i 2 dp, / H + -^(t-t / )l|f|l!-(7;H) 



for almost all t' G J and alH G J with t > t'. 

Using the definition of pj we can rewrite the inequality above as 



-1 f ll u ( t )lL2 + z/ / ll u (s)||wi dsl dpj(u) 



< Pj(U 



I nt')\i,Apj{u) + ^{t-t')ni^Hy 

}(R)) Ju*(R) 



Applying the Monotone Convergence Theorem to the left hand side of the inequality 
above we find, by taking the limit R — > oo, that 



J \\ n ( s )\\h ds } d Pj( u ) 

< f |u(0| 2 L2 dp J (u) + -i-(t-0l|f|| 2 ^(/;H 



for almost all t' G J and allt G J with t > t'. 

Since we are assuming that the mean kinetic energy is finite for almost every i' in 
J, we consider i' which is "good" for the energy inequality and for the corresponding 
finite mean kinetic energy, and such set of t' is still of full measure in /. Then, we 
find from above that 

*-> / Ju(*)|| 2 d Pj (u)GL£ c (J°) 

and is finite everywhere on J, and 

*-> / \\u(t)\\ 2 Hl dp(u) G Ll oc (J°), 
Jul 



where J° is the interior of the interval J. Since these hold for any J C I closed and 
bounded, they imply (3.25) and (3.26). 

Now we go back to the mean energy inequality (3.33) on Uj(kRo) and pass to the 
limit as k — > oo using the the Lebesgue Dominated Convergence Theorem to obtain 
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the mean energy inequality on the whole space Uf. 

jf, {^(luWIi^+^V'duWHOIIu^lll: d S } d Pj (u) 

-L {^(l u (Oli«) + jT^CIuWI^CfW.uW)^ dsj d Pj (u), 
for almost all t' £ J and all t £ J with t >t'. This completes the proof. □ 

Remark 3.4. Notice that the idea of the proof of Theorem 3.1 is to reduce the 
proof of the strengthened mean energy inequality from U\ to the case Uj(R), for 
R > R and J C I compact, and then, for a Vishik-Fursikov measure carried by the 
compact Hausdorff space Uj(R), the idea is to use the Krein-Milman Theorem to 
approximate this measure by a convex combination of Dirac measures concentrated 
on individual weak solutions for which the strengthened mean energy inequality holds. 
This approximation of a Vishik-Fursikov measure by convex combinations of Dirac 
measures will be further explored in Section 4.3 to characterize statistical solutions 
which are projections of Vishik-Fursikov measures. 

Proposition 3.2. Let I C 1R be an interval closed and bounded on the left with 
left end point to. Let p be a Borel probability measure on C\ oc (I,H w ) carried by Uj. 
Suppose that 

[ |u(t)|£ 2 dp(n) < oo, (3.34) 
Ju\ 

for almost all t & I and 

lim / V(|u(t)|| 2 )dp(u)= / ^(|u(t )li 2 )dp(u)<oo, (3.35) 
t-n$ Jul Ju\ 

for every function ip as in Proposition 3.1. Then, 

tH. / |u(t)|i 2 dp(u)£L-(7), (3.36) 

and the strengthened mean energy inequality (3.27) holds including at time t' = t . 

Proof. In view of Proposition 3.1 we only need to prove that the strengthened mean 
energy inequality holds also at t' — to. But this follows trivially from the correspond- 
ing inequality for t' > to and taking the limit as t' — > to using condition (3.35). □ 

Since any Vishik-Fursikov measure satisfies, by definition, the hypotheses in Propo- 
sitions 3.1 and 3.2, we have the following corollary, which we state as a theorem. 

Theorem 3.2. Let p be a Vishik-Fursikov measure over an interval I C K. Then, 
the strengthened mean energy inequality (3.27) holds on I. If I is closed and bounded 
on the left, then (3.27) holds in particular for t' being the left end point of I. □ 
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Remark 3.5. Let / C 1 be an arbitrary interval. For each u G C\ oc (I, H w ) and t G /, 
define 

h + (u,t) = Km sup |u(r)|z,2. 

T-¥t+ 

It is not difficult to see that the strengthened energy inequality (2.8) implies that for 
any uGW, and for any ip as in (2.8), 

^(^ + (u,t)) + ^V / (|u(t)|i 2 )l|u( S )||^ ds 

< ^(/ i 2 + (u,t 1 )) + ^'(|u(t)|i 2 )(f(^),u(s)) L2 ds, (3.37) 
for all t, ti G / with t > t\. Similarly, for a Vishik-Fursikov measure p over J, we find 

|^(^(u,0) + ^V(|u(i)|£ 2 )||u( S )||^ ds} dp(u) 

< {^W(u, + jT^'ClnCOI^Cf uC*))^ da} dp(u), (3.38) 

for all t,ti E J with £ > tj. The advantage of introducing the term h + (u,t) is that, 
with this term, the inequalities hold everywhere, not just almost everywhere and in 
a set depending on the solution. We do not use this inequalities here, but it might 
be useful when considering a collection of weak solutions or statistical solutions. 

3.5. Definition and existence of time-dependent Vishik-Fursikov statistical 
solutions. A Borel probability measure p on C\ oc (I,H w ) induces a family of time- 
dependent Borel probability measures {pt}tei in the phase space H defined by the 
projections p t = HtP, so that 

/ <p(u) dp t (u) = [ <p(v(t)) dp(v), Vt G /, (3.39) 

J H JC loc (I,H w ) 

for every </? which belongs to L l (p t ) for any t G /; in particular for any ip in C^(H W ). 

The following proposition gives the first relation between the statistical solutions of 
the NSE in the sense of Definition 3.2 and the Vishik-Fursikov measures in Definition 
3.3 (see [46, 18]). 

Theorem 3.3. Let I C K be an interval closed and bounded on the left, with the 
left end point denoted to- If p is a Vishik-Fursikov measure over I , then the family 
of projections {pt}t&i is a statistical solution with initial data p t(t (in the sense of 
Definition 3.2). 

Proof. We need to check the conditions in the Definition 3.2 of statistical solution. 
Conditions (ii), (iii), (v), and (vi) follow immediately from the definition of {pt}t&i 
as the projections of p and the properties of p given in Proposition 3.2 and in the 
Definition 3.3 of Vishik-Fursikov. The only remark worth mentioning is that the 
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projection (3.39) from p to p t is initially valid for <p G C b (if w ), but this can be easily 
extended to <f(u) = |u|| 2 , and <p(u) = HuH^, and so on, by using (2.14). 
Condition (V) follows from the fact that 



t ^ / $(u) dp t (u) = / $(u(t)) dp(u) 

Jh Jc loc (i;H w ) 

is the pointwise (in time) limit (thanks to the Lebesgue Dominate Convergence The- 
orem) of the continuous function 



^ / $(P m u(t)) dp(u), 

asm-} oo, where the P m are the Galerkin projectors. 

It remains to prove condition (iv), which is the Liouville-type equation. Let $(u) = 
ip((u, Vi), . . . , (u, Vjt)) be a cylindrical test function as in Definition 3.1. For u G Lij, 
we have t h-> (u(t), Vj) absolutely continuous on I with 

^(u(t),v,) = (F(u(t)),v,V, y G 

where F(u) = f — uAu — B(u, u). Thus, since <fi is continuously differentiable, t t-> 
$(u(t)) is also absolutely continuous on /, with 

d ^ d 

-$(u(t)) = d M^ v 0> • • • > (u, v fc ))-(u(t), V,) 

fc 

= ^(( u > Vl )' • • • ' ( u ' v k))(F(u(t)), Vj) V ',v 

3=1 

= (F(u(t)),&(u(t))) v , >v eL^(I). 

$(u(t)) = $(u(t')) + r(F(u( S )),$'(u( S )))v^,y ds, 

for all t',t G I. Since p is carried by Idj, the previous relation holds p-almost every- 
where, so that, upon integration, 

$(u(t)) dp(u) 

= / $(u(t')) dp(u) + / Af(u(s)), $'(u(s))>v/,v dsdp(u), 

JC loc (I,H v ) </C loc (/,H w ) ./t' 

for all t', tel. The second integrand in the right hand side belongs to L 4//3 (in £ and 
u for the Lebesgue and p measures, respectively) so we apply the Fubini Theorem to 



Hence, 
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write 

*(u(f)) dp(u) 

Cl oc (I,H w ) 

= / $(u(t')) dp(u) + / / (F(u(*)), $'(u( S )))^y dp(u)d S . 

JC loc {I,H v ) Jt' Jc loc (I,H w ) 

Since u i— )■ $(u(t)) is continuous on Ci oc (7, if w ), we have 

/ $(u(t)) dp(u) = / $(u) dp t (u). 

JC loc (I,H v ) JH 

Since the map u !->■ (F(P m u(s)), $'(u(s)))yy is continuous on Ci oc (J, i? w ), and using 
the Lebesgue Dominated Convergence Theorem, we find also that 

f (F(u(s)),<S>'(u(s))} v , y dp(u) 



lim / (F(P m u( S )),$'(u(s))W dp(u) 

lim / (F(P m u),$ / (u))v^, y dp s (u) 
F(u),$'(u)) v ,,v dp s (u). 



Thus, 



/ $(u)dp t (u) = / $(u)dp t ,(u) + / / (F(u),$'(u))v,v dp s (u)ds, 

ifl Jt' JH 

for all t E I. This completes the proof. □ 

There is also the corresponding result for an interval open on the left. 

Theorem 3.4. Let I C R be an interval open on the left. If p is a Vishik-Fursikov 
measure over the interval I, then {pt}tei ^ s a statistical solution on I (in the sense 
of Definition 3.2). 

Proof. Similarly to the proof of Theorem 3.3, except that we need not check condition 
(vi) of Definition 3.2. □ 

A statistical solution in the sense of Definition 3.2 obtained from a Vishik-Fursikov 
measure is called a Vishik-Fursikov statistical solution. More precisely: 

Definition 3.4. Let I C M. be an arbitrary interval. A Vishik-Fursikov statistical 
solution of the Navier-Stokes equations over I is a statistical solution {p t }tei such 
that p t = n t p, for all t £ I, for some Vishik-Fursikov measure p over the interval I. 

The following existence result follows immediately from Theorems 3.1 and 3.3. 
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Theorem 3.5. Let t G E and let po be a Borel probability measure on H satisfying 



Then, there exists a Vishik-Fursikov statistical solution {pt}t>t over the interval I = 



4. Further properties of the Vishik-Fursikov statistical solutions 

4.1. On the carrier of Vishik-Fursikov measures. In this section we are inter- 
ested in the case in which the interval / is closed and bounded on the left. The case 
in which I is open on the left is void of interested for what we are about to discuss 
since in this case the spaces IA\ and Hi are known to be equal. 

We denote the left end point of / by to- We have defined the Vishik-Fursikov 
measures as measures carried by U\ and satisfying in particular a certain strengthened 
continuity at the initial time to of the mean kinetic energy (condition (iii) in Definition 
3.3), although the corresponding energy inequality for the individual solutions in IA\ 
may not be valid at t' = t$. The aim of this section is to prove that a Vishik-Fursikov 
measure is in fact carried by Ui, hence it is carried by the individual solutions for 
which the kinetic energy is continuous at the origin. Notice that since if is a Hilbert 
space and since the kinetic energy is essentially the square of the norm in H, then 
the continuity of the kinetic energy for an individual solution is equivalent to strong 
continuity in H (bearing in mind that weak solutions in Id] are weakly continuous 
also at t = t ). So in essence the result we are about to prove says that the continuity 
of the kinetic energy at time t = to in the mean implies the strong continuity of the 
individual solutions at times t = t almost everywhere. 

Theorem 4.1. Let I C R be an arbitrary interval. Let p be a Vishik-Fursikov measure 
over I . Then p is carried by Ui. 

Proof. As mentioned above, we only need to consider the case in which I is closed and 
bounded on the left, since otherwise both spaces li\ and Uj are known to be equal. 
Consider then such an interval I and let to be the left end point of I. 
Consider the function 



1 - ^0 J to 

which was defined in Lemma 2.5 and proved to be a Borel function on C\ oc (I, H w ) x 
(I \ {t }). Consider also the function 

A(u) = liminf *(u,t), 

which, from Lemma 2.5, is a Borel function in u e «}, with A(u) > for all u e U\, 
and such that A(u) = if and only if u 6 Uj. Therefore, in order to show that p 




[t , oo) satisfying p to = fi . 



□ 




for all t £ I, t > t , 
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is carried by Hi it suffices to show that A(u) = for p- almost all u in U\. Since 
A(u) > for all u G Uj, it suffices to show that 



/ A(u) dp(u) = 0. (4.1) 
Ju\ 



For that purpose, let us recall the mean energy inequality (with ip{r) = r in the 
definition) satisfied by the Vishik-Fursikov measure starting at t' = t^. 



< /^{^|u(to)|! 2 + ^ S (f(r),u(r)) i2 drJ> dp(u) 



for all s G I. Using the Cauchy-Schwarz, Poincare, and Young inequalities on the 
forcing term we find that 

/ t {ns)\h-Hto)\h + v I \Hr)f HX dr\ dp(u)<^( S -* )||f||£~(/ ;J f). 

Taking the time average with respect to s in the interval from to to t and discarding 
the viscous term we find that 

/ *(u,t) dp(u) < ^-(t-^llfll^^), 

for all t G J \ {to}- Let i goes to to to obtain 

limsup / V(u,t) dp(u) < 0. 

Now, since ^(u, t) > — |u(t )|^2 for all t & I and the function u y — |u(t )|^2 is 
p-integrable we may apply Fatou's Lemma to deduce that 

/ A(u) dp(u) = / liminf W(u,t) dp(u) 

Juj(R) Juf(R) t-rto 



< liminf / *(u, t) dp(u) < limsup / *(u, t) dp(u) < 0. 

JuUm t-±t+ Jul 



>U\{R) t->t+ 7w 

Since A(u) > for all u G we deduce from the previous inequality that A(u) = 
for p- almost all u G W|, which means that p is carried by Ui, and the proof is 
complete. □ 
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4.2. Accretion property for Vishik-Fursikov statistical solutions. The evolu- 
tion of the measure of an ensemble of initial data will be studied in the forthcoming 
work [25]. In that work we will be mainly concerned with stationary statistical solu- 
tion, but a particular result about accretion will be first obtained for time- dependent 
Vishik-Fursikov statistical solutions. For that reason, we will mention this result in 
the following remark. 

Remark 4.1. Consider the multi- valued evolution operator defined as follows: Given 
a set E in H and t > 0, we denote by T, t E the set of all points w G H, such that 
w = u(t) and u is in W[o,oo) with initial condition u(0) G E. As it is proved in [25], 
given any Borel set E and any t > 0, the set Yi t E is universally measurable in H, 
i.e. it is measurable with respect to the Lebesgue completion of any Borel measure 
in H. Then, we prove in [25] that if {pt}t>o is a Vishik-Fursikov statistical solution 
over [0, oo), we have that {pt}t>o is accretive with respect to the family {£ 4 } 4 > in 
the sense that 

Pt{V t E) > p (E), 
for all Borel sets E in H and all t > 0. 

4.3. Characterization of the Vishik-Fursikov statistical solutions. After The- 
orems 3.3 and 3.4 one can raise the following question: Given a statistical solution 
{pt}t£i in the sense of Definition 3.2, for some interval 7 C 1, is there a Vishik- 
Fursikov probability measure p on C loc (7~,77 w ) such that U t p = pt, for all t G 77 In 
other words, when can we say that a statistical solution is a Vishik-Fursikov statistical 
solution? We present below a criterion for when this is true or not in the case of an 
interval open on the left and such that pt is carried by a bounded set in H, uniformly 
in t. 

First, let us prove the following localization result. 

Proposition 4.1. A Borel probability measure p onC\ oc (I, H w ) is carried by the space 
C\ oc (I, Bh(R)w) if an d only if pt is carried by Bh{R) for t in a countable dense subset 
of I. 

Proof. Let p be a Borel probability measure on Ci oc (7~, H w ). 

If p is carried by C\ oc (I, Bh{R) w ) and t G / is arbitrary, then, using that IL; 1 B h (R) d 
C\ oc (I , Bh(R)w), we find that 

Pt(B H (R)) = p{Yi^B H {R)) > p(C loc (I,B H (R) w )) = 1, 

so that p t is carried by Bh(R) for all t in /, and, hence, in particular, for t in any 
countable dense subset of I. 

Now let us assume that p t is carried by Bh{R), for t in a countable dense subset 
D of 7". Consider R' > R and let ipR'(s) be the continuous positive function defined 
for s > which is equal to 1, for < s < R, to for s > R', and is linear, for 
R < s < R'. Let <p(u) = ipR'(\P m u\ L 2) for some integer m, where P m is the Galerkin 
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projector defined in Section 2.1. Then, for any time t G D, 

j ?M|P m v(i)| L2 )dp(v)= f tp(y(t)) dp(y) 

<p(u) dp t (u) 



vp R/ (\P m u\ L 2) dp t (u). 
Since p 4 is carried by B H (R) W and ip R >(\P m u\ L 2) = 1 in this ball, we have 
/ V , J R'(l-Pmv(t)| L 2) dp(v) = / ^/j/(|P m u| L 2) dpt(u) 

dp 4 (u) = 1. 



w 

Letting m — > oo, we obtain, thanks to the Lebesgue Dominated Convergence Theo- 
rem, that 

/ ^(|v(t)| ia )dp(v) = l. (4.2) 

JCi oc (I,H w ) 

Letting R' decrease to R, the functions i/jri decrease to the characteristic function 
of the interval [0,-R], so that u i— y V-R'(I u Il 2 ) decreases to the characteristic function 
u \-> xb h (r)( u )- Hence, (4.2) implies, at the limit, 

/ XbhwW*)) dp(v) = l. (4.3) 

Letting 

C/(*;i2) = {v G C loc (I,H w ); \v(t)\# < R} , 
we have that (4.3) implies that 

p{d{t;R)) = 1, Vt G -D. 

Consequently, 



f| =1. (4.4) 



Notice that 

f|C / (t;P) = {vGC loc (J, J PT w ); |v(t)| i2 <i2, VtG £>} . (4.5) 

Now, since t — >■ |v(t)| L 2 is lower-semi-continuous, we have that 
\v(t)\ L 2<R, Vt G R, Vv€ p|C/(t;i2). 
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Thus, 

f)C I (t;R)=C loc (I,B H (R) w ), 

and, hence, 

p(C loc (I,B H (R) w )) = l, (4.6) 
or, in other words, p is carried by C\ oc (I, B H (R) W ). □ 

Now, for any R > Rq, consider the set of statistical solutions (according to Defini- 
tion 3.2) carried by the closed ball Bh{R)'- 

Mr = {{pt}tei is a statistical solution with p t (B H (R)) = 1, Vt G /} . 

One can show that, for such measures, the map 

[ <p(u) d// t (u) eC([t ,oc)), (4.7) 
Jh 

is continuous on /, for any tp G C(Bh{R)w)- Indeed, when p is of the form (3.1), the 
continuous dependence in t of this integral follows from equation (3.2). For a general 
<p G C(Bh(R) w ), we notice that p can be uniformly approximated by functions <pk 
of the type (3.1), thanks to the Stone- Weierstrass Theorem (see Section 3.1). The 
continuous dependence in t of the integrals j H ip>k{xi) d/x 4 (u) gives, at the limit, the 
continuous dependence in t of J H ip k (u) d/x t (u). 

In this set Mr, we define the sequential convergence 

/i*-™-* —>■//, as n — > oo, 

by the condition that, for any ip G C(B H (R) W ), 

p{u) dfi[ n \u) — > / <p(u) d/i 4 (u), as n — > oo, (4.8) 
h Jh 

uniformly in t on any compact subset of /. 

Now, we also define the set of Dirac delta-like statistical solutions, that is statistical 

solutions supported on an individual (Leray-Hopf) weak solution of the Navier-Stokes 

equations, 

— {{$u(t)}tei', u is a weak solution on I with u(t) G B H (R),\/t G /} , 

and the set of Vishik-Fursikov statistical solutions 

VFr = {{p t }t£i is a Vishik-Fursikov statistical sol. with p t (BH(R)) = 1, Vt G /} . 

It is straighforward to show that 

V R c VT R C M r , (4.9) 

and that VFr and Mr are convex. We then consider the convex hull coT>r of T>r, 
which is the set of all finite convex combinations of elements of T>r. Since VFr and 
Mr are convex, we have 

coV R ^VT R ^M R . (4.10) 
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One can also show that VFr and M.r are closed for the sequential convergence 
defined above. Hence, considering the closure of coPr with respect to this sequential 
convergence (i.e. all the measures obtained as the limits of sequences of finite convex 
combinations of measures in T>r), we have 



The next result shows that VFr is precisely the closure of the convex hull of Dr. 
Theorem 4.2. Let I C K be an interval open on the left and R> Rq. Then, 



In other words, a statistical solution in M.r is a Vishik-Fursikov statistical solution 
if and only if it is the limit, in the sense of (4.8), of a sequence of finite convex 
combinations of measures in T>r. 

Proof. It is not difficult to check that T>r is closed in M.r with respect to the conver- 
gence of sequences in A4r according to (4.8). Then, due to (4.9), we need to show 
that 



In other words, we need to show that if a statistical solution belongs to VFr, then 
it can be approximated, in the sense above, by a sequence of convex combinations of 
Dirac measures. Let then {p t }t^i be a Vishik-Fursikov statistical solution such that 
p t = HtP, for all t E I, for some p E C\ oc {I, H w ) with p t (B H (R)) = 1, for all t E I. 

Thanks to Proposition 4.1, we know that p is carried by C\ oc (I, Bjj(R) w ), and, 
hence, it is carried by U\{R). 

Since U\{R) is a compact Polish space, and thus a compact separable space, and p 
is a Borel probability measure on Uf(R), it follows from the Krein-Milman Theorem 
that each p is the limit of a sequence of convex combinations of Dirac deltas in Uj(R), 
i.e. 




(4.11) 



VF R = (coV r ). 




(4.12) 




(4.13) 



with of > o, J2j = X > v f e u i(R)- 



The convergence (4.13) means that, for every $ E C(U}(R)), we have 




j 



Taking $(v) = <p(v(t)), for <p E C(B H (R)) and t E I, define 



Wf) = £*i n V(vj B) ) 
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and notice that 



Mt) = E *j B VM B) ) /" v(v(*)) d P (v 



V?(u) d/>t(u) = ^(t), Vt G /. 



Since the weak solutions in Ll\{R) are weakly equicontinuous on any compact in- 
terval J G I with values in H w , we have that the if) m are equicontinuous, hence their 
pointwise convergence to ip on J implies their uniform convergence to ip on J, for 
every compact interval J C I. This means that {X^^j v-^ft)}* conver g es to {p*}* 
in the sense of (4.8), which concludes the proof of (4.12) and of Theorem 4.2. □ 

We now have the following characterization of the Vishik-Fursikov statistical solu- 
tions. 

Corollary 4.1. Let I C M be an interval open on the left. Let \i = {fi t }tei be 
a statistical solution in Mr, for some R > Rq. Then, {fi t }tei is a Vishik-Fursikov 
statistical solution if and only if there exists a sequence {p^jn of convex combinations 
of statistical solutions of the form {5 v (t)}tG/; where v e U\{R), such that — > \i in 

Mr. 

Proof. Since any family of the form {8v(t)}tei * s a Vishik-Fursikov statistical solution, 
any convex combination of them, as in the statement of the corollary, is a Vishik- 
Fursikov statistical solution as well, and, therefore, so is its limit. This last fact as 
well as the converse property follow directly from Theorem 4.2. □ 

Remark 4.2. In the case of measures with unbounded support, not belonging to 
Mr for any R > 0, we consider, for E > 0, the set 

M(e) = |{/-it}te-T is a statistical solution ; J |u|| 2 dp t (u) < E, Vt G l\ . (4.14) 

We also define in a similar manner the sets "Die) an d VjF<e)- Then, 

V (E) C coV {E) C VF^) CM {E ), 

and the analog of Theorem 4.2 would be that every statistical solution in M( E ) is 
a Vishik-Fursikov statistical solution if and only if it is the limit, in some suitably 
defined sense, of measures in co T)<e)- However, this analog result remains a conjecture 
at this point and depends, in part, on an appropriate definition of convergence for 
such measures. □ 
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